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ABSTRACT 

We have incorporated the latest release of the Padova models into the evolutionary synthesis 
code Starburst99. The Padova tracks were extended to include the full asymptotic giant branch 
(AGB) evolution until the final thermal pulse over the mass range 0.9 to 5 M©. With this ad- 
dition, Starburst99 accounts for all stellar phases that contribute to the integrated light of a 
stellar population with arbitrary age from the extreme ultraviolet to the near-infrared. AGB 
stars are important for ages between 0.1 and 2 Gyr, with their contribution increasing at longer 
wavelengths. We investigate similarities and differences between the model predictions by the 
Geneva and the Padova tracks. The differences are particularly pronounced at ages > 1 Gyr, 
when incompleteness sets in for the Geneva models. We also perform detailed comparisons with 
the predictions of other major synthesis codes and found excellent agreement. Our synthesized 
optical colors are compared to observations of old, intermediate-age, and young populations. 
Excellent agreement is found for the old globular cluster system of NGC 5128 and for old and 
intermediate-age clusters in NGC 4038/39. In contrast, the models fail for red supergiant dom- 
inated populations with sub-solar abundances. This failure can be traced back to incorrect red 
supergiant parameters in the stellar evolutionary tracks. Our models and the synthesis code are 
publicly available as version 5.0 of Starburst99 at http://www.stsci.edu/science/starburst99/. 
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Starburst99 is a multi-purpose evolutionary 
synthesis code for modeling various properties of 
stellar populations. The code was first used by 
Leitherer & Heckman (1995) for compiling a set 
of observable parameters of evolving young stellar 
populations. The Leitherer & Heckman models 
relied on an earlier release of stellar evolutionary 
tracks by Maeder (1990). These tracks were su- 
perseded by a new set of models from the Geneva 
group (Meynet et al. 1994 and references therein) 
so that an update to the original synthesis code 
seemed appropriate. This was done by Leitherer et 
al. (1999), who implemented the latest set of the 
Geneva tracks, added the capability of perform- 
ing isochrone synthesis, and released a web-based 



interface to offer users interactive model computa- 
tions atwww.stsci.edu/science/starburst99/. Sev- 
eral important additional aspects were addressed 
over the years: Leitherer, Robert, & Drissen 
(1992) discussed the mechanical energy release 
of evolving stellar populations; Leitherer, Robert, 
& Heckman (1995) studied the satellite-ultraviolet 
(UV) spectra at high spectral resolution; de Mello 
et al. (1998) extended the latter work to less mas- 
sive B stars; Leitherer et al. (2001) added an em- 
pirical UV library of metal-poor stars; Robert et 
al. (2002) included 0-star spectra obtained with 
the FUSE satellite; most recently, Smith, Norris, 
& Crowther (2002) implemented the latest gener- 
ation of blanketed non-LTE model atmospheres. 

Starburst99, like any other evolutionary synthe- 
sis code, can in principle be applied to stellar pop- 
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ulations having arbitrary properties. In practice, 
the useful parameter space that can be rehably 
predicted is determined by the set of stellar at- 
mospheres, spectral templates, and evolutionary 
tracks. Starburst99 relies on the most advanced 
atmospheres and spectra for hot, massive stars 
(Lejeune, Cuisinier, & Busser 1997; 1998; Hillier 
& Miller 1998; Schmutz 1998; Pauldrach, Hoff- 
mann, & Lennon 2001) and their coupling to the 
Geneva evolutionary models (Schmutz, Leitherer, 
& Gruenwald 1992). This property makes Star- 
burst99 and similar codes (e.g., Schaerer & Vacca 
1998, Cerviiio, Mass-Hesse, & Kunth 2002 and 
references therein) uniquely qualified for model- 
ing populations containing hot massive stars, in 
particular Wolf-Rayet (W-R) stars. On the down- 
side, Starburst99 is not competitive when it comes 
to stellar populations whose properties are dom- 
inated by intermediate-age and old stars. This 
regime is better described by the stellar evolu- 
tion models of the Padova group (Fagotto et al. 
1994b and references therein) whose input physics 
pays particular attention to intermediate- and low- 
mass stars. Furthermore, the Padova models offer 
the possibility to extend the post-main-sequence 
evolution of intermediate-mass stars until and in- 
cluding thermally pulsing asymptotic giant branch 
(TIP-AGB) stars. These properties make the 
Padova tracks the models of choice for synthesis 
codes studying populations with ages above ^ 10^ 
years. The most important synthesis codes con- 
centrating on this age range are, e.g., those of 
Worthey (1994), Vazdekis et al. (1996), Fioc & 
Rocca-Volmerange (1997), and Bruzual & Char- 
lot (2003). 

Over the years, it has become customary to 
model young and old stellar populations with the 
codes that were optimized for the parameter space 
of interest and perform a comparison between the 
codes in the overlap region at intermediate age 
where either model is supposed to be applicable 
(e.g. Kewleyetal. 2001; Whitmore & Zhang 2002; 
Lamers et al. 2002). Very often, however, such 
a comparison is complicated by the fact that ad- 
ditional ingredients (such as model atmospheres) 
vary between different codes. Therefore our group 
has expanded the capabilities of Starburst99 to al- 
low for a self-consistent treatment of stellar phases 
starting with and past the early asymptotic giant 
branch (EAGB) phase. In this paper we describe 



the upgrade done in Starburst99 and how it affects 
the model predictions. Since Starburst99 now of- 
fers models based on both the Geneva and the 
Padova tracks, relative comparisons can be done 
in a reliable way. 

An additional major goal of this paper is a crit- 
ical evaluation of the reliability of our evolution- 
ary synthesis models when red supergiants (RSG) 
dominate. Stellar evolution models are notori- 
ously uncertain for RSGs, in particular at low 
metallicity (Massey & Olsen 2003). Therefore, 
comparisons of synthesis models with extragalac- 
tic systems containing massive cool stars are badly 
needed to provide empirical guidance for stellar 
evolution models. 

This paper is organized as follows. First, we re- 
view similarities and differences between the main 
astrophysical ingredients in both the Geneva and 
Padova evolution models (Section 2). In Section 3 
we investigate how those ingredients affect the 
main stellar predictables from both sets in the 
Hertzsprung-Russell diagram (HRD). A technical 
discussion of implementation issues in Starburst99 
is in Section 4. Our main focus during the imple- 
mentation was the treatment of the AGB stars. 
We discuss our approach in Section 5. A compre- 
hensive comparison between the predicted prop- 
erties of a standard population using either the 
Geneva or the Padova tracks in Starburst99 is per- 
formed in Section 6. In Section 7 we extend this 
comparison to other synthesis codes. Comparisons 
of old, intermediate-age, and young populations 
are performed in Section 8. Finally, in Section 9 
we provide the conclusions, as well as guidelines 
for using the new models in Starburst99. 

2. Main physical ingredients in the Padova 
and Geneva tracks 

In this section we will compare the major phys- 
ical ingredients in the two sets of evolutionary 
models. This assessment will help us identify the 
root causes of the differences in the predicted stel- 
lar properties of the two model sets, and eventu- 
ally in the differences of the population proper- 
ties synthesized by Starburst99. Throughout this 
paper we will use the "high mass loss" tracks of 
Geneva, as opposed to the "standard" tracks. The 
former are the recommended set whose properties 
have been fine-tuned to match the observed prop- 
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erties of massive stars (Mcynct ct al. 1994; but sec 
Massey 2003 for an opposing point of view). The 
relevant literature sources for the Geneva mod- 
els arc Schallcr ct al. (1992), Schacrcr ct al. 
(1993a,b), Charbonnel et al. (1993), and Meynet 
et al. (1994). 

The Padova models used here are those for 
intermediate and low-mass stars in the range of 
0.15-7 M0 from Girardi et al. (2000) and the 
tracks for massive stars in the range of 9 - 120 Mq 
from the original set of Padova models discussed 
in Bressan et al (1993), Fagotto et al. (1994a,b), 
and Girardi et al (2000). Throughout this work we 
will refer to this composite set of tracks as "Padova 
models" . 

We will focus our discussion in this section 

on the adopted chemical composition, the nuclear 
opacities and reaction rates, the convective energy 
transport, the mass-loss rates, and the boundary 
conditions at the stellar surface. 

2.1. Chemical composition 

Both sets of tracks use the solar-abundance ta- 
bles of Anders & Grevesse (1989) as a basehne, 
with updated values for carbon, nitrogen, and iron 
from Grevesse & Anders (1991), Grevesse ct al. 
(1991), Hannaford et al. (1992), and Grevesse & 
Noels (1993). Note that these abundances are still 
on the older solar oxygen abundance scale with 
log(0/H) + 12 = 8.9. More recent work by Al- 
lende Prieto, Lambert, & Asplund (2002) has led 
to a downward revision of the solar oxygen abun- 
dance. The currently preferred value is log(0/H) 
+ 12 = 8.69. 

Both the Padova and the Geneva models adopt 
the Anders & Grevesse values for the zero- age- 
main-sequence (ZAMS) abundances of all ele- 
ments other than He. (We will discuss the he- 
lium abundances further below.) In Table 1 we 
summarize these values for the most important el- 
ements up to iron. In the following we will adopt 
these conventions and nomenclature: abundances 
of individual elements, when the elements are de- 
noted by their chemical symbols, are always by 
number on a scale where the hydrogen abundance 
equals 12. Following the nomenclature used in 
stellar evolution, we will use X, Y, and Z for the 
mass fractions of hydrogen, helium, and all ele- 
ments heavier than helium, respectively. We will 



avoid the incorrect use of "mctallicity" to denote 
Z. Metallicity should be restricted to the logarith- 
mic number abundance of metals relative to the 
sun. 

The heavy-element composition for non-solar 
abundance models in both the Geneva and the 
Padova tracks was derived from the solar-abundance 
models by applying an appropriate scaling of the 
heavy-element abundances and by keeping the 
metal-abundance ratios fixed. Small adjustments 
to the metal-abundance ratios were required in or- 
der to be consistent with the heavy-element ratios 
used for the opacity calculations (Schaller et al. 
1992). Both Padova and Geneva published five 
complete sets of evolutionary tracks with chem- 
ical compositions covering a wide range of as- 
trophysical applications. Three of the five sets 
have identical metal abundances in Geneva and 
Padova: Z = 0.02 (solar), 0.008 (40% solar), and 
0.004 (20% solar). In addition, Padova pubhshed 
super-solar models with Z = 0.05 (2.5 solar) for 
the full mass range and with Z = 0.03 (1.5 so- 
lar; Girardi et al. 2000) with updated opacities 
for intermediate- and low-mass stars. Our com- 
posite set of Padova tracks uses a merged ver- 
sion of the Z = 0.05 and 0.03 tracks for massive 
and intermediate-|-low-mass stars, respectively. 
For comparison, Geneva published models with 
Z = 0.04 (twice solar). Padova and Geneva 
also published extremely metal-poor tracks with 
Z = 0.0004 (2% solar) and 0.001 (5% solar), re- 
spectively. The differences in the heavy-element 
abundances should be kept in mind when making 
comparisons between different tracks. 

The helium abundances follow from the solar 
and primordial helium abundances, respectively, 
and the adopted AY/AZ relation. The Padova 
models use a primordial helium abundance of 0.24 
(note that Girardi ct al. 2000 use 0.23) and 
AY/AZ = 2.5. The Geneva group prefers a 
slightly steeper AY/AZ = 3 for sub-solar and so- 
lar models, whereas for super-solar models they 
use AY/AZ = 2.5, in agreement with Padova. 
Although the helium abundance is important in 
controlling the opacity and determining the core 
temperature profile, the effects of varying Y be- 
come dominant only at high metal abundance (Yi 
2003). 

A compilation of the ZAMS abundances of H, 
He, C, N, and O for the five metal abundances is in 
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Table 2. In summary, we find very similar chem- 
ical abundances for both the Geneva and Padova 
models. Except for the very metal-rich and -poor 
model sets, different abundance profiles cannot be 
the reason for differences between the tracks. 

2.2. Nuclear opacities and reaction rates 

The opacities in both sets of evolutionary mod- 
els are from the widely used Livermore group 
(e.g., Iglesias, Rogers, & Wilson 1992 and ref- 
erences therein). The reaction rates adopted by 
the Padova and Geneva models are different. The 
largest disagreement is for the rate of ^^C(a, 7)^^0 
in the CN cycle. The Geneva group developed its 
models with the rates calculated by Caughlan et 
al. (1985), whereas the Padova group used those 
of Caughlan & Fowler (1988). The latter calcula- 
tions result in a lower reaction rate. 

The effects of varying the rate of this reaction 
on evolutionary models have been known since the 
early study by Iben (1972). The larger cross sec- 
tion leads to a higher conversion rate of carbon 
into oxygen and to more pronounced loops ex- 
tending to the blue part of the HRD before rapid 
core contraction and envelope expansion set in, 
and the evolution proceeds back to the red. Over- 
all, the core-He-burning lifetime is increased by 
several percent. In practice, however, these effects 
are difficult to separate from those of mass loss 
and convective overshoot (Chiosi 1998). 

2.3. Convective energy transport 

A comprehensive discussion of the importance 
of the convective energy transport for stellar evolu- 
tion models and synthetic stellar populations was 
given by Yi (2003). While Yi's work is tailored to- 
wards intermediate- and low-mass stars, many of 
the arguments are also qualitatively correct when 
applied to massive stars. In particular, the effi- 
ciency of the convective energy transport has pro- 
found consequences on the evolution of massive 
stars. Massive stars have convective cores and ra- 
diative envelopes, except for the W-R phase. In- 
creasing the size of the convective core, e.g., by 
convective overshooting, will increase the rate of 
the energy transport since convection is more ef- 
fective than radiation. Therefore the surface tem- 
perature will critically depend on these assump- 
tions. 



The prescription for the convective overshoot- 
ing in the Padova tracks is as follows: low- mass 
stars have a core overshoot of d/H-p = 0.25, where 
d is the overshoot distance and Hp is the pressure 
scale height, whereas more massive stars have a 
value of 0.5. In contrast, this parameter is set 
to 0.2 in the Geneva models, as suggested by an 
empirical calibration over the mass range from 
1.25 to 25 Mq (Schaller et al. 1992). Fagotto 
et al. (1994a) demonstrated that the Padova and 
Geneva overshoot parameters are in fact consis- 
tent, as the algorithm for calculating the overshoot 
distance is different in the two sets of evolution 
models. 

The mixing length parameter a = I /Hp, with 
I being the mixing length, was derived by the 
Geneva group from an empirical calibration of 
their models via the observed location of the red 
giant branch in 75 star clusters (Schaller et al. 
1992). a = 1.6 ±0.1 was found and adopted. The 
Padova models use a very similar mixing length 
parameter of o; = 1.63 (Bressan et al. 1993). 

Both sets of models make rather similar as- 
sumptions on the convective energy transport. 
However, neither set considers the effects of rota- 
tion. Rotation has been recognized to be a major 
driver of stellar evolution, in particular for mas- 
sive stars (Meynet & Maeder 2003, and references 
therein). Systematic models suitable for imple- 
mentation into synthesis models are under con- 
struction. When such models will become avail- 
able, their post-main-sequence evolution may dif- 
fer from that in the current Padova and Geneva 
models. 

2.4. Mass-loss rates 

Stellar mass loss is a fundamental parameter for 
the evolution of massive stars. Both Geneva and 
Padova use the parameterization of de Jager et 
al. (1988) who compiled empirical mass-loss rates 
across the entire upper HRD. Deviating from this 
prescription, the Geneva group modified the pre- 
scription for the mass-loss rates of W-R stars and 
their evolutionary progenitors. All massive stars 
with radiatively driven winds (i.e., OB and WNL 
stars) are assumed to have twice the de Jager 
mass-loss rates. The motivation for this devia- 
tion from de Jager's relation is purely empirical: 
the resulting W-R parameters in the HRD are in 
better agreement with the observations. The ad- 
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justment of the mass-loss rate is simply a proxy 
for other effects not accounted for, such as ro- 
tation and convection. The rates of WNE and 
WC stars were calculated with the prescription 
of Langer (1989), who showed that these stars 
tightly follow the mass-loss rate vs. mass rela- 
tion M = (0.6-1.00) X lO^(Af/A^0)2-5 [Mq yr-^]. 
For the Padova tracks, the W-R phase sets in for a 
hydrogen surface abundance below 30% (Bressan 
et al. 1993), whereas for the Geneva tracks this 
limit is 40% (Schaller et al. 1992), in both cases 
the mass-loss rate from Langer's relation is used. 

The mass-loss rates of stars with non-solar 
chemical composition were scaled with the relation 
M oc ^0-5 (Kudritzki, Pauldrach, fc Puis 1987). 
This scaling was not applied in the W-R phase 
for the Geneva models. In contrast, the Padova 
tracks include this scaling of the mass-loss rates 
in all stellar phases. 

Neither set of evolutionary models includes 
mass loss for low- and intermediate-mass stars, 
as the rates are too low to be significant. When 
we added the thermally pulsing AGB stars, we 
adopted the Reimers relation (Rcimcrs 1975) for 
all evolved stars in the red part of the HRD with 
masses below 12 Mq outside the AGB phase. 

2.5. Boundary conditions at the stellar 
surface 

The most significant difi^erence between the 
Padova and the Geneva tracks occurs in the W-R 
phase, when the definition of effective temperature 
(Teff) becomes ambiguous because of the highly 
extended atmospheres. The challenge is to define 
the radius and Tcff in the atmosphere. The calcu- 
lation depends on the choice for the opacity in the 
wind of the W-R star. Langer (1989) and Maeder 
(1990) proposed to calculate the opacity from the 
continuity equation assuming a standard velocity 
law V{r) oc (1 — R/rY , where R is the radius of 
the star and r is the distance. The simplifying as- 
sumption of pure electron scattering opacity (cTg) 
is made. The radius refers to the location where 
the optical depth t{Rcs) = 2/3, and Toff is deter- 
mined from the Stefan-Boltzmann law at -Rcff- 

In the case of the Geneva tracks, Schaller et 
al. (1992) followed Langer's (1989) suggestion, to- 
gether with the Castor, Abbott, and Klein the- 
ory (Castor et al. 1975; Pauldrach et al. 1986; 



Kudritzki et al. 1989) to calculate the opacity. 
They adopted a mean weighted opacity k = ae{l + 
FM) obtained from the radiation-driven wind the- 
ory, where FM is the so-called force- multiplier, 
which gives the contribution from absorption lines. 
This parameter is calculated from the radiation- 
hydrodynamics in the atmosphere. The final i?cfr 
and Teff are calculated by adding the line opacity 
to the pure electron scattering opacity. This cor- 
rection for Toff shifts the stellar positions to the 
red part of the HRD. Effective temperatures ob- 
tained with this correction cover the domain of the 
observed WNL, WNE, and WC stars shown as by 
Maeder (1990) and Schaerer & Maeder (1992). 

In the Padova tracks, there is no correction 
for the temperature decrease due to stellar winds. 
The boundary conditions simply follow from the 
optical depth calculated with the opacities of 
Rogers & Iglcsias (1992). Therefore, no outflow 
of material is allowed for, and all stars follow the 
same definition of Toff. As a result, W-R stars in 
the Padova models have higher temperatures than 
those in the Geneva models. 

3. Physical differences and their effects on 
the stellar evolution 

In this section we will analyze the effects of 
the different physical input parameters on the 
main evolutionary properties of individual stars 
in the HRD. First, we will discuss the princi- 
pal parameters luminosity L and Teff as a func- 
tion of time. Then we will describe the differ- 
ences in the surface abundances. We have cho- 
sen those three sets of tracks whose abundances 
are identical in the Padova and Geneva models 
{Z = 0.004,0.008,0.02). Both sets of tracks pre- 
dict very similar stellar parameters (except in the 
W-R phase), with the solar metallicity tracks from 
the Padova group producing somewhat older ages 
than those of the Geneva group (~ 10% and ~ 4% 
for intermediate and massive stars, respectively). 

On the ZAMS, L and Toff are very similar for 
most stellar masses in both model sets. Padova 
models are slightly hotter by ~ 0.03 dex and more 
luminous by ^ 0.02 dex than Geneva models. The 
similarity of the values is due to the same source 
of the opacities and abundances, and the similar 
treatment of overshooting. As discussed by Char- 
lot, Worthey, & Bressan (1996) and Yi (2003), the 
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major effect is the overshooting, yet the opaci- 
ties and the element abundances play an impor- 
tant role as well. Bruzual & Chariot (2003) noted 
the 50-200 K higher temperature of the red gi- 
ant branch (RGB) in the models of Girardi et al. 
(2000) in comparison with the original Padova set. 
Consequently, the Girardi et al. tracks produce 
bluer colors than those from the original Padova 
set in the Bruzual & Chariot models after the ap- 
propriate color transformation is applied. 

The evolution of L and T^s as a function of time 
is very similar, except for the W-R phase where 
the models from Geneva are more luminous and 
cooler than the models from Padova. This behav- 
ior is expected due to the different boundary con- 
ditions on the stellar surface in the Geneva models. 
The ultimate astrophysical reason for the differ- 
ence is not so much a real difference in T^e but 
rather a different definition of temperature. T^s 
in the Padova models refers to deeper atmospheric 
layers, whereas the Geneva definition samples the 
outer atmospheric region where the photons can 
escape (Schmutz, Lcithcrcr, & Grucnwald 1992). 

In the case of intermediate-mass stars, we note 
that both sets of models do not take into ac- 
count the tip of the asymptotic giant branch (TIP- 
AGB). In order to minimize confusion, we avoid 
calling the end of the AGB phase the first thermal 
pulse in the AGB phase (TP- AGB) as opposed to 
the convention often used in other work. The be- 
havior of Tcff is very similar in Geneva and Padova 
except when the horizontal branch (HB) is mod- 
eled. The Padova tracks have a different treat- 
ment for this phase. As a result, the temperature 
becomes higher than in the Geneva tracks. The 
behavior of these parameters in the AGB phase 
will change when we introduce the TIP- AGB for 
the Padova tracks as discussed in Section 5. 

For the low-mass regime, the set of tracks from 
the Geneva group does not include the EAGB 
phase for stars with 0.8 < M/A1q < 1.7, whereas 
the Padova tracks include the evolution up to the 
EAGB phase for stars M > 0.6 Mq. StarburstQQ 
does not use the supplemental release of low-mass 
tracks by the Geneva group (Charbonnel et al. 
1996). This release extends the early evolution of 
stars in the mass range 0.8 < M/Mq < 1.7 to the 
EAGB phase. As a result, the Starburst99 im- 
plementation of the Geneva models in this mass 
range lacks complete coverage of even the RGB 



phase. Figure 1 shows tracks for the low-mass 
regime emphasizing the different phases covered 
by each group. The 1.7 Mq track is the transi- 
tion between stars with and without a degenerate 
helium core in the Geneva models and are repre- 
sented with a solid line and crosses, respectively, 
in Fig. 1. The track without a degenerate core has 
a clump of points which represent the counterpart 
of the EAGB phase in the Geneva tracks. This 
transition has been discussed by Girardi & Bertelli 
(1998). The additional evolutionary phase in the 
extended 1.7 Mq Geneva track (and the tracks 
corresponding to lower masses) will produce no- 
ticeable features in the evolution of colors. 

We also compared the evolution of He, C, N, 
and O as a function of time. For massive stars, 
the values reached at the beginning of and during 
the W-R phase are very similar at the three abun- 
dances. The final value of the oxygen abundance 
is ~ 0.3 dex higher in the Geneva tracks than in 
the Padova tracks (e.g., at 60 and 40 Mq). The 
carbon abundance is very similar in both groups, 
with a difference of ~ 0.1 dex at 60 and 40 Mq. 
For the same masses, the nitrogen abundance be- 
haves similarly, except for models at Z = 0.008 
where the Padova 40 Mq model has a much higher 
value than the Geneva model. This is the result 
of different limiting masses for the formation of 
W-R stars. These differences do not significantly 
affect the integrated parameters of synthetic pop- 
ulations. However, they do affect the census of 
W-R subtypes, which are defined by the surface 
abundance of these elements. 

For intermediate stellar masses, wc found abun- 
dance differences of < 0.01 dex in boths sets of 
tracks. 

4. Implementation of the Padova tracks in 

Starburst99 

We implemented the publicly available set of 
Padova models into Starburst99, closely following 
the scheme wc used for the Geneva tracks. The 
models are distributed as tables with the perti- 
nent stellar parameter given as a function of time 
for each stellar mass. The main differences be- 
tween the Padova and the Geneva models, in addi- 
tion to the previously identified W-R parameters, 
are the larger mass range and the more detailed 
coverage of late evolutionary phases of low- and 
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intermediate-mass stars in the Padova tracks. 

Two sets of Padova models were prepared for 
implementation. One set has the stellar evolu- 
tionary phases as specified in the Padova tables, 
and the second includes additional entries for ther- 
mally pulsing AGB stars. Apart from these addi- 
tional entries, the two sets are identical, and both 
are available to users of Starburst99. We will dis- 
cuss the specifics of the AGB phase in the subse- 
quent section. 

As wc did in the Geneva tracks, wc selected 
the initial mass Mj, stellar age T, current mass 
Mc, logi, logTeff, and the surface abundances of 
H, He, C, N, and O by mass fraction from the 
Padova tables. The original tables list M only for 
massive stars with masses > 9 Mq. The missing 
entries were compiled as follows. For the AGB 
phase, we used the values of Vassiliadis & Wood 
(1994). For stars off the main-sequence with miss- 
ing table entries, we applied the Reimcrs (1975) 
formula. Finally, we adopted low place holder 
values of 10~^^ M© yr~^ for less massive main- 
sequence stars. This choice is purely driven by 
computational reasons. Main-sequence stars with 
masses below 9 Mq are unimportant for the mass 
return of a population and, in addition, their mass- 
loss rates are poorly known. We performed these 
updates of the mass-loss rates only in the set of 
Padova tracks to which we added the TIP- AGB 
stars. The original Padova release was left un- 
changed. Note the absence of the core temper- 
atures and radii, which are used in the Geneva 
tracks to distinguish wind-corrected W-R atmo- 
spheres from their uncorrected counterparts. We 
set the core temperatures and radii equal to Tcfr 
and i?eff, respectively, in our implementation of 
the Padova tracks. Furthermore, we assumed the 
same mass limits for the formation of W-R stars 
in the Padova models as in the Geneva models. 
The lowest stellar mass evolving into W-R stars 
depends on the metal abundance, as suggested by 
observations (Massey 2003). The Geneva values 
are Mi = 21, 25, 35, 42, and 61 Mq for twice so- 
lar, solar, 40%, 20%, and 5% solar composition, re- 
spectively (Maedcr & Meynet 1994). We adopted 
the same values for the Padova tracks with chem- 
ical composition 2.5 solar, solar, 40%, 20%, and 
2% solar, respectively. 

A successful interpolation technique requires 
identification of the strategic evolutionary phases 



along the tracks. These are: the ZAMS point, 
the first temperature minimum prior to core- 
hydrogen exhaustion, the terminal-age main- 
sequence (TAMS), the epoch of core-hydrogen 
exhaustion, the supernova phase, the beginning 
of the RGB phase and the epoch of its luminosity 
maximum in the HRD, the onset of the horizontal 
branch evolution and its maximum luminosity, the 
epoch of core-helium exhaustion, and the begin- 
ning of the EAGB phase. During the interpola- 
tion, special attention is paid to these phases to 
ensure no unphysical result is obtained for those 
initial masses which are not explicitly specified in 
the distributed Padova tables. 

After linking the Padova tables to Starburst99, 
they are treated identically to the Geneva tracks 
(see Leitherer et al. 1999 for details on the com- 
putational procedure). Either a traditional mass 
interpolation or an isochronc interpolation in time 
is performed. The latter technique is particular 
suitable for later evolutionary phases which would 
otherwise produce discontinuities in the synthe- 
sized quantities. Extensive test were performed to 
ensure the interpolation produces sensible results. 
We are confident that the implementation of the 
Padova models into Starburst99 with and without 
the AGB phases is as reliable as that of the Geneva 
tracks. 

In addition to allowing for the AGB phases, 
the Padova tracks extend the valid mass range in 
Starburst99 from the previous 0.8 Mq limit in the 
Geneva tracks down to 0.15 Mq in the Padova 
tracks. Not only does this remove the previous re- 
striction on modeling low-mass stars in young pop- 
ulations. More significantly, it allows Starburst99 
to model stellar populations of any age up to the 
oldest cosmological ages with the best available 
stellar models. 

5. AGB models: first thermal pulse up to 
the planetary nebula 

We extended the published Padova tracks be- 
tween 0.9 and 5 Mq to include thermally pulsing 
AGB stars. In this section we describe the adopted 
stellar properties in the AGB phase and their im- 
plementation in the Padova tracks. 

The evolution of the AGB phase after the first 
thermal pulse has been the subject of numer- 
ous studies (e.g., Groenewegen & de Jong 1993; 
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Blocker 1995; Dorfi & Hofner 1998; Ventura 2001; 
Wachter et al. 2002). Parameters of interest 
are the high mass-loss rate, variability in tem- 
perature and luminosity, and the time scale of 
this evolutionary phase. Observational guidance 
is provided, e.g., by Mira variables, where vari- 
ability has been observed and attributed to ther- 
mal pulses (Wood & Zarro 1981). The mass loss 
in this stage has been related to planetary neb- 
ula (PN) nuclei as well as to white dwarfs (WD) 
whose mass distribution peaks around 0.6 Mq. 
Therefore, WDs are considered the evolved de- 
scendants of AGB stars (Wcidemann & Koester 
1983; Weidemann 1990). The progenitors of AGB 
stars are assumed to be > 1.0 Mq (Pottasch 1983; 
Jura 1990). The upper mass limit for AGB for- 
mation is somewhat uncertain but lies between 
5 and 8 Mq, the latter value representing the 
transition to the formation of a core-collapse su- 
pernova. Renzini & Voli (1981) noted that the 
mean mass-loss rate required to produce a typical 
PN was 3 X 10"^ Mq yr~^ They introduced 
the term superwind to describe such high mass- 
loss rates since they are much higher than those 
of RSGs described with the formula of Reimers 
(1975). Since then, many AGB stars (particularly 
infrared sources, OH/IR stars, and Mira variables) 
have been found to have mass-loss rates consistent 
with the superwind expectation (e.g., Knapp & 
Morris 1985 in the Galaxy, and Wood et al. 1992 
in the LMC). 

Given the lower mass limit of the progenitor of 
~ 1.0 Mq, it is clear that AGB stars are important 
for galactic ecology. The amount of gas mass con- 
verted into newly formed stars in this mass range 
is ~ 25% of the total, depending on the initial 
mass function (IMF). Therefore, these stars be- 
come important for the total stellar mass loss in 
systems with a significant intermediate-age popu- 
lation. 

The central issue of the relation between the 
mass loss during the AGB evolution and other 
parameters, including the core mass of the star 
has received wide attention (e.g., Groenewegen & 
de Jong 1993; Vassihadis & Wood 1993, Blocker 
1995). The formulation of Groenwegen & de Jong 
has been applied to the Padova tracks before (e.g., 
Girardi & Bertelh 1998; Marigo 1998, Girardi et 
al. 2000). Since the Reimers (1975) mass-loss re- 
lation fails to reproduce AGB winds, an empirical 



approach relating M and other periodic proper- 
ties in the thermally pulsing AGB phase was de- 
veloped by Vassihadis fc Wood (1993). We chose 
their models to develop a set of TIP- AGB parame- 
ters and add them to the evolutionary tracks from 
Padova. 

Vassihadis & Wood (1993) published models 
with Z = 0.004, 0.008, and 0.02 for the mass range 
0.9 - 5.0 M0. Their models follow the evolution 
from the first thermal pulse in the EAGB phase 
up to the stripping of the core in the PN (or TIP- 
AGB) and WD phases. Their definition of the end 
of the AGB phase is the time when Toff becomes 
higher by 0.3 dex than the value for a "reference" 
AGB star prior to blueward evolution. In our im- 
plementation of the Vassihadis & Wood models we 
terminated the evolution of the AGB phase in the 
final pulse before the last ejection of material at 
the lowest value of Teff. 

The models give Mc, Tcs, L, and M during 
the evolution from the first AGB pulse until the 
last pulse prior to the ejection of the nebula. We 
have interpolated these tracks in the range (0.9 - 
5.0 M0) to be consistent with the Padova tracks. 
The tracks were also re-sampled to cover the same 
evolution phases with a small number of points. 
The Padova tracks were re-sampled from the point 
where = up to the TIP- AGB phase (the pulse 
immediately prior to the lowest temperature) with 
17 points over the tracks including the first ther- 
mal pulse already modeled in tracks. We used the 
Reimers formula to calculate the mass loss during 
the preceding RGB phase. 

A complication arises from the lack of AGB 
models for Z = 0.0004 and 0.05. It is currently 
unknown if and how the chemical composition af- 
fects the AGB evolution. In the absence of other 
evidence, we assumed the Z = 0.004 AGB mod- 
els to be applicable to the Z = 0.0004 Padova 
tracks, and the Z = 0.02 models to the Z = 0.05 
tracks. No further adjustment to the AGB models 
was made prior to the implementation. The main 
properties of the final set of tracks are summarized 
in Table 3. In this table we contrast the imple- 
mentation of the Geneva and the Padova tracks 
with AGB stars in Starburst99. We list the range 
of abundances, mass range, the number of grid 
points, and the main evolutionary phases. Se- 
lected tracks with different chemical composition 
are shown in Fig. 2. In this figure we have high- 
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lighted the extension of the original Padova tracks 
towards cooler temperatures and higher luminosi- 
ties with the addition of the TIP-AGB phase. As 
we will demonstrate below, this phase has signif- 
icant effects on the predicted properties of stellar 
populations. In Fig. 3 we have plotted three tracks 
showing the AGB mass range at higher resolu- 
tion for different chemical composition. Included 
in this figure are the Geneva tracks, the original 
Padova tracks, and the Padova tracks with the 
TIP-AGB phase. The comparison highlights the 
differences between the three sets of tracks for in- 
termediate and low stellar masses. Lower metal 
abundances shifts all tracks to higher tempera- 
tures. In addition, the TIP-AGB correlates in- 
versely with L: over most of the AGB mass range, 
more metal-poor stars are more luminous than 
their metal-rich cousins. 

6. Starburst99: Geneva vs. Padova tracks 

In this section we provide a comparison of syn- 
thetic populations generated with Starburst99 in 
order to identify similarities and differences intro- 
duced by the new set of tracks. A comparison 
with other synthesis codes will be performed in 
Section 7. 

The model simulations in this and the following 
section were performed (unless otherwise stated) 
for a standard single stellar population of mass 
10^ M© with a Salpeter IMF between 0.1 and 
100 Mq. We chose to use a low- mass cut-off of 
0.1 Mq for the IMF, instead of the canonical value 
of 1 M© . While this choice introduces an offset by 
a factor of ~ 2.5 with respect to the original Star- 
burst99 models, it allows us to include low-mass 
stars below 1 M© in the simulations, which would 
otherwise be excluded. 

We used the fully line-blanketed atmospheres of 
Smith, Norris, & Crowther (2002) for hot stars and 
ATLAS9-|-Phoenix models as adjusted by Lejeune, 
Cuisinier, & Buser (1997; 1998) for stars of spec- 
tral type B and later. The ATLAS9+Phoenix at- 
mospheres were also used to assign spectral energy 
distributions (SED) to stars in the AGB phase. 
While this simplification is adequate for modeling 
SEDs, dedicated AGB spectra, such as those of 
Mouhcine & Langon (2002), are required for the 
synthesis of spectral lines. 

The stellar atmospheres as well as most auxil- 



iary tables in Starburst99 have chemical compo- 
sitions tailored for the Geneva models. Since the 
Padova and Geneva tracks have different chemi- 
cal composition at the low {Z ~ 0.001 and 0.0004 
for Geneva and Padova, respectively) and the high 
[Z = 0.04 and 0.05, respectively) extremes, there 
is a slight mismatch in Z. We used the avail- 
able data files (i.e., with Z = 0.001 and 0.04) 
for the most metal-poor and -rich Padova tracks. 
Consequently, the synthesized properties for these 
tracks have larger intrinsic errors than those for 
the Z = 0.004, 0.008, and 0.02 tracks where our 
abundances arc consistent. 

Previously published Starburst99 models were 
all performed with linear time steps. While this 
technique works fine for young ages, logarithmic 
time steps are more appropriate for populations 
with ages of order Gyr. Otherwise, either the time 
resolution at young ages would be unacceptably 
low, or the data volume would be unacceptably 
high. Consequently, the release of Starburst99 
containing the Padova models offers the choice of 
linear or logarithmic time steps. In the following 
we will discuss representative model simulations 
with both the Geneva and Padova tracks for sin- 
gle populations. We also calculated corresponding 
models with continuous star formation but will not 
plot the results because of redundancy with the in- 
stantaneous models. 

6.1. Spectral energy distributions 

We have compared the stellar SEDs for the full 
age range. Throughout this paper we follow the 
convention to plot the Geneva, original Padova, 
and Padova+TIP-AGB models with solid, dashed, 
and dotted lines, respectively. Young populations 
with solar composition observed at 1, 3, 10, 50, 
and 100 Myr arc shown in Fig. 4. SEDs produced 
with the Padova tracks are more luminous by up 
to 0.1 dex at most wavelengths. These differences 
are the result of higher intrinsic stellar luminosi- 
ties in the Padova models for post-main-sequence 
evolutionary phases. 

There is a striking difference between the SEDs 
in the UV around 3 Myr. The excess emission in 
the Padova models is produced by hot W-R stars 
when coupled with extended model atmospheres. 
While this coupling is consistent when done for the 
Geneva models, it produces spurious results in the 
case of the Padova models. When W-R stars are 
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important, the preferred approach is to use the 
Geneva tracks together with the extended atmo- 
spheres. If the Padova tracks are used, a minimal- 
istic technique would be to couple them with the 
static ATLAS9 atmospheres. We will discuss the 
ionizing continuum in more detail in Section 6.4. 

We selected five different ages of 0.1, 0.5, 1, 
5, and 12 Gyr as representative for intermediate- 
age and old stellar populations. These models 
are plotted in Fig. 5. The contribution of the 
TIP-AGB stars to the red part of the spectrum 
is clearly visible: for models with ages > 1 Gyr 
(e.g., 5 Gyr), there is a strong difference between 
the Padova and Geneva tracks. The difference re- 
sults from a combination of two effects. First, the 
Geneva tracks do not include low-mass stars be- 
low 0.8 Mq. Since our chosen IMF extends down 
to 0.1 Mq, the absence of these stars creates a lu- 
minosity deficit, which sets in at an age of a few 
Gyr, and becomes progressively stronger. Second, 
the implementation of the Geneva models in Star- 
burst99 does not include the add-on tracks pub- 
lished by Charbonnel et al. (1996). Charbonnel 
et al. computed the horizontal branch evolution 
after the He-flash in low-mass stars. Horizontal 
branch stars around ~ 1.5 Mq contribute signifi- 
cantly to the near-IR luminosity beginning at an 
age of ~ 1 Gyr. They are the chief reason for 
the difference between Geneva and Padova at and 
after that epoch. 

6.2. Massive-star inventory 

The number of O stars derived from Padova 
models is the same as that from the Geneva mod- 
els. This confirms the consistency between the two 
sets of tracks for the O phase in massive stars. The 
slight luminosity and temperature difference in the 
tracks is too minor to affect the stellar census in 
Starburst99. 

Similarly, the supernova rates following from ei- 
ther set of tracks is essentially the same since the 
stellar lifetimes are very similar and the cut-off 
mass for supernova formation is the same. 

In Fig. 6 we have plotted the W-R/ O star ratios 
predicted by the Geneva and the Padova tracks. 
The ratios agree quite well, as expected from the 
similar W-R definitions in both sets. Recall that 
W-R stars are defined via their surface abundances 
and (for computational reasons) via their lower 



mass limit for formation. We apply the same defi- 
nitions for both sets of tracks in Starburst99. Ob- 
viously, had we tied the W-R definition to T^s, the 
agreement in Fig. 6 would largely disappear be- 
cause of the different temperature definitions. The 
small differences between the W-R/ O ratios in Fig. 
6 are caused by slightly different stellar lifetimes 
and numerical noise. As a reminder, W-R/0 is 
strongly dependent on the chemical composition. 
A decrease in metals inhibits the formation of W- 
R stars. 

W-R stars come in two subclasses: nitrogen- 
rich WN and carbon-rich WC stars. Their classi- 
fication is tied to the surface abundances of He, C, 
N, and O. Details are in Leitherer et al. (1999). 
The WC/WN ratio is an important IMF tracer 
(Massey 2003). Fig. 7 suggests reasonable agree- 
ment between the predictions resulting from the 
Padova and the Geneva tracks. On average, the 
Geneva models lead to higher WC/WN ratios. 

6.3. Ionizing continuum 

We compared the ionizing fluxes from hot, mas- 
sive OB stars resulting with the new sets of tracks. 
No significant differences are found for the hydro- 
gen ionizing photons at any of the flve chemical 
compositions. The fluxes shortward and close to 
the Lyman edge are predominantly contributed by 
OB stars, whose number densities are very similar 
for the two sets. W-R stars are negligible contrib- 
utors for a standard population mix. 

W-R stars do matter for the ionizing fluxes at 
shorter wavelengths in the He° and He"*" contin- 
uum below 504 and 228 A, respectively. The differ- 
ent Teff definition in the Padova and Geneva tracks 
causes enormous differences in phases when W-R 
stars are present. This is illustrated in Figs. 8a and 
b for the number of photons capable of ionizing 
He° and He+, respectively. In order to highlight 
how different model atmospheres interact with dif- 
ferent tracks, we present four atmosphere/track 
combinations in the two figures. "Atm 2" refers 
to the plane-parallel Lejeune atmospheres in Star- 
burst99, whereas "Atm 5" is the atmosphere pack- 
age released by Smith et al. (2002), which as- 
signs fully blanketed models accounting for stellar 
winds. 

Prior to the appearance of the first W-R stars 
(at ~ 3 Myr), the four atmosphere/track combi- 
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nations predict Hc*^ ionizing photons which agree 
within a factor of 2 (Fig. 8a). During the W- 
R phase, the predictions diverge by up to a fac- 
tor of 10. Wc consider the combination of the 
Geneva tracks and Smith et al. (2002) atmo- 
spheres as the currently best available models for 
W-R related parameters and use them as our base- 
line (solid line in Fig. 8a). When the Geneva 
tracks are coupled to the Smith atmospheres in 
Starburst99, wc do not use the surface tempera- 
ture but a hotter value from deeper in the atmo- 
sphere. This hotter temperature is the required 
input for the spherical extended atmospheres. For 
comparison, linking the Lejeune atmospheres with 
the Geneva tracks leads to much lower ionizing 
fluxes from W-R stars because we use the lower 
surface temperatures in the evolution models and 
because the atmospheric fluxes of given tempera- 
ture are slightly different from those of Smith ct 
al. (long-dashed line). Contrasting these results 
with those from the Padova tracks, one recognizes 
that the higher Tog in the Padova models always 
leads to higher ionizing fluxes, irrespective of the 
atmosphere. The small diff'erences between the 
results for the static (dash-dotted) and expand- 
ing (short-dashed) atmospheres coupled with the 
Padova tracks are again caused by the diff'erent 
atmospheric fluxes. 

The interplay between the tracks and the at- 
mospheres becomes even more dramatic at wave- 
lengths below 228 A (Fig. 8b). O stars dominate 
the luminosity until about 3 Myr. The Smith 
et al. (2002) atmospheres use the WMBASIC 
wind models of Pauldrach et al. (2001) for O- 
star phases. These wind models predict 3 to 4 or- 
ders of magnitude higher luminosities than static 
ATLAS9 (Lcjcunc) models for both the Geneva 
and Padova tracks (solid and short-dashed lines 
in Fig. 8b). When W-R stars become important. 
Smith et al. rely on the CMFGEN atmospheres 
of Hillier & Miller (1998). At solar and higher 
composition, the wind models predict recombina- 
tion of He+ in the W-R envelopes, and the emer- 
gent fluxes below 228 A are quite low (solid line). 
Since the wind recombination is not accounted for 
in the static ATLAS9 atmospheres, the prediction 
from the "Geneva-I-Atm 2" model (long-dashed) 
becomes unrealistically large. The behavior of the 
two Padova models is as expected: the dominant 
effect is the higher T^g, which boosts the He+ 



luminosity by orders of magnitude. Because of 
the very high temperatures, recombination is less 
important, and the extended atmospheres predict 
only somewhat lower luminosities than the static 
models with the Padova tracks (short-dashed and 
dash- dotted lines). 

The Padova models show a pronounced deflcit 
of radiation around ages of 5 Myr when linked to 
the extended atmospheres. This depression results 
from the very short lifetime of a 25 M© W-R star 
in the models. The 25 Mq track is unique in that it 
has a hot W-R star only during the last few thou- 
sand years prior to the termination of the track. 
In contrast, the 40 Mq model predicts the hot 
W-R phase to last for several hundred thousand 
years. Since the extended W-R atmospheres are 
only linked to the evolutionary models when both 
a minimum temperature and a minimum helium 
surface abundance are exceeded, there will be a 
period of a few 10^ yr around 5 Myr (correspond- 
ing to ZAMS masses around 30 Mq) without any 
stars that would satisfy the definition required by 
the atmospheres. 

The graphs in Figs. 8a and b become very dif- 
ferent at lower and higher chemical composition. 
This results from the combination of the depen- 
dence of the W-R formation on chemistry and the 
behavior of the W-R atmospheres on mass loss and 
metal abundance. Smith et al. (2002) give an ex- 
tensive discussion of these effects as predicted for 
the Geneva models. As the Padova models are not 
primarily designed for modeling W-R stars, ex- 
tending Fig. 8 to other chemical abundances would 
not add useful information to the discussion in the 
present paper. 

6.4. Panchromatic luminosities 

We compare the predicted luminosities at rep- 
resentative wavelengths in Fig. 9. As before, the 
Geneva, original Padova, and Padova-|-TIP-AGB 
models are presented. Fig. 9 shows the bolomet- 
ric (Mboi), 1500 A (Li5oo), 5 band (Mb), Vhand 
(My), J band (Mj), and /('band (Mk) luminosi- 
ties. The passband definitions are the same as 
in Leitherer et al. (1999). All filters are in the 
Johnson (1966) system except for R and /, which 
are in the Cousins system (Bessell 1983; 1990). 
-^1500 is the average of the luminosity between 
1490 and 1510 A, and MboI is referenced to a sys- 
tem in which the sun has MboI = 4.75. (The lat- 
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ter zero point is equivalent to using the relation 
log(L/i0) = -0.4Mbo1 + 1.90.) 

The overall trend in Fig. 9 suggests very good 
agreement betwcicn the three models until about 
100 Myr when the first AGB stars form. The AGB 
phase lasts from 100 Myr to 2 Gyr and produces 
excess luminosity mostly in the IR. At epochs later 
than about 3 Gyr, the Geneva and Padova models 
diverge at all wavelengths, including MboI. This 
is caused by the same cfFccts already discussed in 
Section 6.1: the Geneva tracks do not account for 
the post-He-flash evolution in intermediate-mass 
stars and they lack low-mass stars altogether. As 
a result, the Geneva tracks lead to a light deficit 
when these stars become important. 

The readers and users of Starburst99 should be 
aware of uncertainties inherent in L1500 for ages 
older than ~ 3 Gyr. After that epoch, UV-faint 
F and later-type stars provide the UV luminos- 
ity. The standard Kurucz model atmospheres are 
a poor match for the UV light of old populations 
(Lotz, Ferguson, & Bohhn 2000). Efforts are un- 
derway to remedy the situation (Peterson, Dor- 
man, & Root 2001). 

The results for different chemical abundances 
arc qualitatively similar to those shown in Fig. 9 
for solar composition. Changing the composition 
has its biggest effect in the RSG dominated phase 
when the evolutionary tracks predict strong varia- 
tions of the RSG temperatures and lifetimes with 
abundance. We will test this prediction in Sec- 
tion 8.3. 

The fractional contribution of AGB stars to 
the luminosity in four passbands is reproduced 
in Fig. 10. Plotted is the magnitude increase by 
adding AGB stars relative to the prediction by the 
original Padova tracks. The four panels in Fig. 10 
show MboI, Mb, My, and Mk for three differ- 
ent compositions. As suggested by the trends in 
Fig. 9, AGB stars have minor influence on optical 
passbands, such as B and V, but can contribute 
up to 1 mag to near-IR filters. The contribution 
to the bolometric luminosity is 0.4 mag or less. 
This behavior is rather sensitive to variations of 
chemical composition as a result of the models of 
Vassiliadis and Wood (1993). 



6.5. Colors 

The color evolution is of particular observa- 
tional interest, as it is a normalized quantity and 
does not depend on the (a priori unknown) star- 
formation history. We have plotted our results for 
{B - V), {V - J), and {V - K) in Fig. 11. In this 
figure, we include the color evolution for models 
with Z = 0.02, 0.008, and 0.004 in order to identify 
abundance effects. The overall trends are similar 
to those in Fig. 9. The Geneva and Padova models 
agree, except during periods when stellar phases 
dominate which are not well accounted for or are 
not included at all in one model set or the other. 
This is the case during the first few Myr when 
W-R stars contribute (not relevant to Fig. 11), 
during the AGB phase, and after ~ 1 Gyr when 
post-Hc-flash stars become important. The RSG 
dominated phase around 10 - 20 Myr differs in 
the Geneva and Padova tracks because of differ- 
ent lifetimes and temperatures. We will return to 
this issue in Section 8. 

As a rule of thumb, main-sequence stars are al- 
ways important contributors to optical and near- 
IR colors, except during the RSG dominated pe- 
riod between 10 and 20 Myr (the large depres- 
sion in Fig. 11). AGB stars become compara- 
ble in their luminosity around 1 Gyr. At even 
older ages, RSGs provide the bulk of the luminos- 
ity. These arguments are of course dependent on 
the actual color and should be viewed as a rough 
guideline only. See also Chariot & Bruzual (1991) 
and Maraston (1998) for a related discussion. 

Fig. 11 nicely illustrates how the colors are sen- 
sitive to chemical composition: more metal-poor 
stars are generally bluer than metal-rich stars (ex- 
cept during the AGB phase - see below). This 
is the joint result of bluer positions in the HRD 
and a steeper decrease of line- blanketing at shorter 
wavelengths for lower abundances. The RSG de- 
pression at ^ 10 Myr has a striking Z dependence 
in both sets of tracks. There are fewer and warmer 
RSGs at lower chemical abundances, a theoreti- 
cal prediction discussed by Cervifio & Mas-Hesse 
(1994). We will perform a series of observational 
tests of this prediction in Section 8. 

The colors of AGB dominated phases are sensi- 
tive to chemical composition as well. The effect is 

most pronounced in the red and IR. This compo- 
sition dependence is introduced by the Vassiliadis 
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and Wood models, which predict a higher AGB lu- 
minosity at lower abundance. Consequently, AGB 
stars (which are very red) contribute more light at 
low Z. There is mild evidence for this effect from 
a comparison of Galactic and Magellanic Cloud 
AGB stars (Vassihadis & Wood 1993). 

6.6. Mass loss and mass-to-light ratio 

Starburst99 calculates quantities related to the 
mass and energy return in addition to the photon 
output. In the following we will use the mass- loss 
related properties as they follow from the table 
entries in the evolutionary tracks in order to be 
able to compare the tracks themselves. This differs 
from the recipe in Leitherer et al. (1999), where 
the stellar-wind properties were derived from a hy- 
brid theoretical and empirical formula, indepen- 
dently of the stellar evolution models. 

The mass loss by winds from all newly formed 
stars is plotted in Fig. 12. The predictions fol- 
lowing from either the Padova or the Geneva 
models are essentially identical until the appear- 
ance of the first AGB stars. The difference be- 
tween the Geneva and Padova models around 
100 Myr is caused by the lack of M entries for 
stars with M < 12 Mq in the Padova tables. 
The Padova-|-AGB tracks extend the prediction 
for M beyond ~ 1 Gyr, at which epoch Star- 
burst99 would otherwise cease to be useful. AGB 
stars with their strong winds are the prime source 
for the mass return to the ISM for intermediate- 
mass stars. Nevertheless, the total mass lost in 
the AGB phase is minuscule when compared to 
the mass loss in earlier stellar phases. This can 
be appreciated in Fig. 13, where we have plotted 
the cumulative mass loss for the same stellar pop- 
ulation. The AGB contribution is the incremental 
mass loss increase noticeable after about 3 Gyr. 

The corresponding mass-to-light ratios are in 
Fig. 14. We show the predictions for LboI, Ly, and 
Lk- Other passbands behave accordingly. The 
three sets of tracks lead to very consistent results, 
except for those phases with missing stellar ingre- 
dients, as discussed before. We emphasize that 
the results in Fig. 14 do not take into account the 
removal of stellar material during the evolution 
of the population. In other words, M is constant 
and taken to be the initial total stellar mass. In 
practice, M decreases because some (or even most) 
of the mass lost by stellar winds can escape from 



the gravitational field of the stellar population and 
become unbound. Therefore, the mass-to-light ra- 
tios in Fig. 14 are lower limits to the astrophysical 
situation. 

7. Comparison of Starburst99 with other 
evolutionary synthesis codes 

How do the predictions of Starburst99 for the 
photometric properties of stellar populations with 
and without AGB stars compare to those from 
other synthesis codes? We surveyed the literature 
and collected the relevant data from those sources 
which either published the color evolution directly 
or allowed us the extract the photometry from the 
information provided in the publication. 

The literature sources used in this comparison 
are listed in Table 4, together with their principal 
ingredients. Included in this table are only those 
sources which specifically address the AGB evo- 
lution. Buzzoni (1989) released a widely used set 
of synthetic population properties but his mod- 
els commence after 1 Gyr. Therefore we did not 
include them in our comparison. Col. 1 of this ta- 
ble gives the model designation. The evolutionary 
models are in col. 2. The source for the AGB evo- 
lution is in col. 3. In col. 4 we give the upper mass 
limit Mmax included in the AGB evolution. The 
adopted SEDs and the photometric system are in 
cols. 5 and 6, respectively. In col. 7 we provide 
the literature source. Our point of reference is the 
Starburst99 implementation of the Padova-|-TIP- 
AGB models (this work). The second and third 
rows in Table 4 list the recent "GALAXEV" re- 
lease by Bruzual & Chariot (2003) who updated 
their earlier work by linking the new intermediate- 
resolution empirical library of Le Borgne et al. 
(2003). The GALAXEV models include the TIP- 
AGB evolution for both the original and the Gi- 
rardi versions of the Padova tracks. Except for the 
use of the empirical library, the ingredients in the 
GALAXEV models are closest to those of Star- 
burst99 among the entries in Table 4. "Spectral" 
(Vazquez, Carigi, & Gonzalez 2003) incorporates 
both the Padova isochrones as released in the pa- 
per of Girardi et al. (2000) or the Padova (1994) 
isochrones (Bertelli et al. 1994), respectively^. 

The Padova isochrones should not be confused with the 
tracks. The isochrones contain T^g and L, and these quan- 
tities are directly used by Spectral for the synthesis 
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The Spectral ingredients include a very approxi- 
mate model for TIP-AGB stars from Boothroyd 
& Sackmann (1988) and Groenewegen & de Jong 
(1993) in the 1994 and 2000 isochrones, respec- 
tively. Mouhcine & Langon (2002) published ded- 
icated AGE evolution models, based on the ear- 
lier work of Wagenhuber & Groenewegen (1998). 
Maraston's (1998) models are based on the fuel- 
consumption theorem (Buzzoni 1989) and incor- 
porate the stellar evolution models of Castellani, 
Cheffi, & Straniero (1992). These evolution mod- 
els are for the assumption of no overshooting, mak- 
ing them a useful test for this particular ingredi- 
ent. For the AGB phase the models from Blocker 
and Schonberner were used. Finally, we include 
the "PEGASE" code of Fioc & Rocca-Volmerange 
(1997) which is complete for the full AGB phase, 
including the post- AGB evolution as calculated by 
Schonberner (1993) and Blocker (1995). 

In addition to differences in the adopted evo- 
lution models, the eight codes in Table 4 consider 
quite different mass ranges for the AGB phase. We 
note that the lower boundary mass is only of aca- 
demic interest as long as it is below 1 Mq because 
less massive stars have not had time yet to reach 
the AGB phase since the first stars formed in the 
universe. The upper boundary mass is more rele- 
vant. Maraston (1998) assumes the highest value, 
i.e., 8 Mq, whereas most other work uses 5 Mq. 
The Starburst99 value is 5 Mq. Fortunately, these 
differences have little effect on the colors of the 
population. Massive AGB stars, while luminous, 
are quite rare because of the IMF weighting, and 
they have shorter life times. Therefore the rele- 
vant masses are between 3 and 5 M©, which are 
covered by all entries in Table 4. 

The synthesis codes in Table 4 fall into two 
groups with respect to their assumptions for 
the SEDs. Fioc & Rocca-Volmerange (1997), 
Mouhcine & Langon (2002), and Bruzual & Char- 
lot (2003) incorporate empirical library stars to 
reproduce AGB stars, the other models (includ- 
ing Starburst99) rely on stellar atmospheres. The 
SED entries in Table 4 are high-level guidelines 
only. Often the authors use an elaborate mix 
of several libraries with cross-calibrations. The 
salient point is that one group of entries in Ta- 
ble 4 is primarily based on empirical, and the 
other on theoretical SEDs. Finally, the photomet- 
ric systems can differ. Some predictions are in the 



original Johnson (1966) system, whereas others 
are in the Johnson-Cousins system as defined by 
Bessell (1983; 1990) and Bessell & Brett (1988). 
While the different filter profiles and zero points 
are small in most filters, they do matter for the R 
and / passbands. 

None of the listed codes uses exactly the same 
ingredients as Starburst99. Nevertheless, by com- 
paring the different model predictions on a rela- 
tive basis, we can identify the trends introduced 
by different assumptions and attempt to verify the 
consistency of Starburst99 with other work. 

In Fig. 15 we have plotted the color evolu- 
tion in the {U - B), (B - V), {V - R), and 
{y — K) colors for solar abundance. For the sake 
of clarity we plotted Starburst99 and GALAXEV 
in the left panels, and the remaining models in 
the right panels. (Note that Mouhcine & Langon 
2002 did not publish optical colors. Their IR 
results are discussed separately further below.) 
Except for random numerical noise, Starburst99 
and GALAXEV have no systematic differences by 
more than ^ 0.05 mag in the optical colors, and 
~ 0.2 mag in the near-IR. We do not detect signif- 
icant differences with respect to either the Padova 
(1994) or the Girardi (2000) implementations in 
GALAXEV. Although this agreement was hoped 
for because of the similar inputs, it is nevertheless 
a gratifying result. 

The comparison with the predictions of Spec- 
tral, Maraston (1998), PEGASE, and Mouhcine & 
Lancon (2002) is in the right panels of Fig. 15. We 
note the generally very good agreement in [U — B) 
and {B — V) at any age and in iy—R) and iV—K) 
outside the AGB phase. When AGB stars dom- 
inate, the predictions differ by up to one mag- 
nitude. The differences can be understood from 
the earlier discussion of the different stellar tracks 
and AGB models used by the different codes. The 
predictions of Maraston (1998) in particular differ 
becaiisc of their use of a fully independent set of 
stellar evolution models. The agreement between 
the predictions of the other codes does not nec- 
essarily suggest those models arc preferable; the 
codes responsible for those predictions are simply 
using the same evolutionary tracks. Furthermore, 
Maraston's {V — R) colors are in the original John- 
son (1966) system. The R filter has a pivot wave- 
length of 0.69 /um, whereas the Cousins R is at 
0.64 ^,m. This leads to redder {V — R) colors in the 
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Maraston models. The Padova isochroncs in their 
implementation in Spectral lead to very high AGB 
luminosities. This is the result of the adopted 
temperature-luminosity relation hard-coded in the 
Padova isochrones release. 

Next, we turn to the near-IR colors. Since 
some literature sources provide only the near-IR 
colors, others only the optical colors, and some 
a mix of both wavelength ranges, we are pre- 
senting the near-IR results separately in Fig. 16. 
This figure serves a dual purpose. First, we are 
comparing the AGB sensitive near-IR colors for 
different codes, and second, we are evaluating 
the influence of the heavy-element abundance. 
The comparison at solar chemical composition 
is consistent with the previous results at shorter 
wavelengths. Starburst99 and GALAXEV are in 
excellent agreement. The comparison with the 
models of Vazquez ct al. (2003), Fioc & Rocca- 
Volmerange (1997), and Mouhcine & Langon 
(2002) is somewhat less favorable but neverthe- 
less still acceptable. The trend with composition 
in the AGB phase is opposite to that seen at older 
ages: lower abundances lead to redder colors. 

Finally, we provide a comparison of the mass- 
to-light ratio predicted by Starburst99 using the 
Padova-l-TIP-AGB tracks and by those sources in 
Table 4 which have published this quantity. The 
mass-to-light ratios in Fig. 17 were obtained un- 
der two assumptions for the mass evolution in 
Starburst99. In one case we adopted a constant 
mass equal to the initial total stellar mass, and 
in the other we allowed for a continuous mass de- 
crease according to the mass loss prescription in 
the stellar-evolution models. We are distinguish- 
ing between these two cases for a fair comparison 
with other models, some of which were computed 
only with (Maraston 1998) and some only without 
(Fioc & Rocca-Volmerange 1997) mass loss. The 
GALAXEV models consider both cases. Fig. 17 
suggests excellent agreement overall. Starburst99 
and GALAXEV are the best matches for either 
mass-loss model. As expected, Starburst99 is in 
better agreement with Fioc & Rocca-Volmerange 
when M is constant in Starburst99 (upper panels). 
On the other hand, the models of Maraston agree 
better with Starburst99 for non-conservative evo- 
lution (lower panels). Again, we note the impact 
of the different stellar evolution models used by 
Maraston and Starburst99. 



8. Testing the predictions of Starburst99 
for intermediate-age and old popula- 
tions 

In this section we will perform basic perfor- 
mance tests of our models by comparing the pre- 
dicted colors to sets of relevant observational data. 
As before, we will ignore very young, ionizing pop- 
ulations since the new set of evolutionary models is 
not optimized for these populations, and extensive 
tests of Starburst99 applied to hot stellar popula- 
tions have been reported in the literature. Three 
age regimes will be explored: ^ 10^° yr, as traced 
by old globular clusters, ~ 10* yr observed in su- 
per star clusters, and ~ 10^ yr, the age regime of 
RSGs. 

8.1. Globular clusters in NGC 5128 

Peng, Ford, & Freeman (2004a, b) performed an 
in-depth study of the globular cluster system sur- 
rounding NGC 5128. It is not our goal to repeat 
their excellent analysis. Rather, we will take ad- 
vantage of their photometric database of 215 con- 
firmed globular clusters with radial velocities and 
compare the measured {U — B), (B — V), {V — R), 
and {V — I) colors to those predicted by Star- 
burst99. This comparison is done in Fig. 18 for 
{B — V) and (V — I). Wc constructed two-color 
diagrams for the other colors as well, with consis- 
tent results. These diagrams are not reproduced 
here. 

The theoretical colors in Fig. 18 cover the ago 
range from to 10 Gyr, with tick marks high- 
lighting 0.1, 1, 5, and 10 Gyr. Models for three 
chemical abundances are plotted: Z — 0.0004, 
0.004, and 0.02. We restrict our discussion to 
models using the Padova-|-TIP-AGB tracks since 
the Geneva tracks are known to fail at those old 
ages. The models illustrate a textbook example of 
the notorious age-metallicity degeneracy (Worthey 
1994): the metallicity and the age vectors are al- 
most perfectly parallel at old ages, making it im- 
possible to separate the two parameters in the ab- 
sence of additional information. The data points 
are the entries from the catalog of Peng et al. 
(2004a). They are not corrected for foreground 
reddening. Peng et al. (2004b) favor a redden- 
ing of E{B — V) = 0.115, and the corresponding 
shift of the data points is indicated by the redden- 
ing vector. Note that the vector is parallel to the 
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direction of the metallicity and age evolution. 

The colors predicted by Starburst99 coincide 
with the {B — V)-{V — I) plane defined by the ob- 
servations. As discussed by Peng et al. (2004b), 
the globular cluster system of NGC 5128 consists 
of two well-defined age and metallicity groups, one 
having {B - V)o = 0.71 and {V - I),, = 0.89, and 
the other 0.88 and 1.10, respectively. (The sub- 
scripts indicate dereddened colors.) The mean col- 
ors of the two cluster groups arc in excellent agree- 
ment with, e.g., the age > 10 Gyr for Z = 0.004 
and the age < 10 Gyr for Z = 0.02 models. Of 
course other age/metallicity combinations in the 
models reproduce the data as well. The main point 
is that the models are an excellent match to the 
data, given the restrictions imposed by the age- 
metallicity degeneracy. 

8.2. Super star clusters in NGC 4038/39 

NGC 4038/39 ("The Antennae") arc the archety- 
pal interacting galaxies whose newly formed 
system of super star cluster is interpreted as 
an important piece of support in favor of the 
merger-formation hypothesis for elliptical galax- 
ies (Schweizer 1998). Whitmore et al. (1999) ob- 
tained deep HST WFPC2 photometry of NGC 4038/39 
and identified distinct age groups ranging from 
very young ionizing clusters, over intermediate-age 
10* yr), to old globular clusters. We can uti- 
lize their published photometry for a model com- 
parison with an emphasis on intermediate ages. 
We proceed as in the case of NGC 5128, i.e., we 
constructed two-color diagrams for all published 
passbands but we will reproduce only the {V — I) 
vs. (B — V) diagram (Fig. 19). The lay-out of 
Fig. 19 is identical to that of Fig. 18. We show the 
same sets of models, and we include the published 
data points without the appropriate reddening 
correction. The latter is rather substantial for 
young clusters because of the significant internal 
reddening of E{B - V) = 0.3 in NGC 4038/39 
(Whitmore et al. 1999). 

Whitmore et al. (1999) detected 11 very red 
clusters which they interpreted as old globular 
clusters from the original pre-merger galaxies. 
Comparison with the models in Fig. 19 confirms 
their suggestion. The 11 clusters occupy the same 
parameter space as the globular clusters of NGC 
5128 in Fig. 18. Any of the three plotted age 
sequences would be able to reproduce their colors 



with a suitably chosen age/metallicity combina- 
tion. 

The intermediate-age clusters fall into the age 
range 0.25 to 1 Gyr according to Whitmore et al. 
(1999). 25 such clusters were identified by Whit- 
more et al. and have been added to Fig. 19. Al- 
most all super star clusters in this age group arc 
clustered aroimd (B-V) ~ 0.25 and (V-I) ~ 0.4. 
This is the age when the first AGB stars influence 
the c;oIors in our models. The models are not af- 
fected by an age-metallicity degeneracy in this age 
range (~ 200 Myr) because of the dependence of 
the AGB phase on chemical composition. It is im- 
portant to understand the physical reasons for the 
vertical loops at (S - F) ~ 0.0 to 0.3 in Fig. 19. 
The loop for the Zq models is caused by RSGs, 
which are much less important lower abundances. 
The redder {V — I) at higher Z reflects the lower 
Tcff of RSGs for higher abundances. In contrast, 
the loops in the Z = 0.004 and 0.0004 models are 
caused by AGB stars, which are predicted to have 
higher L at lower Z. Since both RSGs and AGB 
stars are very red, they are much more important 
for red than for blue stars. Therefore, {B — V) 
does not show these Z dependences, and the dis- 
tinct loops form in Fig. 19. NGC 4038/39 has 
solar chemical composition, and the correspond- 
ing model sequence is a rather good match to the 
data. Conversely, the models with sub-solar chem- 
ical composition are clearly excluded. Of course, 
this prediction calls for independent verification 
with more data, given the results with find for 
RSGs in Section 8.3. Nevertheless, we can state 
that models with appropriate abundance do in fact 
reproduce the observed cluster colors. 

The youngest clusters in the sample of Whit- 
more et al. (1999) are characterized by the bluest 
colors and their association with HII regions. 
These 50 objects in Fig. 19 exhibit the largest 
dispersion. The most likc;ly cause is strong con- 
tamination by nebular and [O III]. The V 
photometry was obtained by transforming the 
original WFPC2 F555W data into the Johnson 
system. The F555W filter encompasses several 
strong nebular lines, including the mentioned 
and [O III]. In contrast, the F439W (used in lieu 
of B) and the F814W (used in lieu of T) filters 
do not suffer from such contamination. Therefore 
the data points representing the youngest clus- 
ters in Fig. 19 must be corrected towards bluer 
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{B — V) and redder {V — I), i.e., towards the up- 
per left. In addition, the clusters are affected by 
strong internal reddening, which is not accounted 
for. A quantitative interpretation of the youngest 
clusters is not our goal. Instead we will test the 
model predictions for stellar populations with ages 
of tens of Myr with blue compact dwarf galaxies 
observed in the IR. 

8.3. The red stellar population of blue 

compact dwarf galaxies 

Our next test will specifically address the va- 
lidity of our evolutionary synthesis models for 
younger age ranges when RSGs arc important. 
This is roughly the epoch between 8 and 30 Myr 
(see Fig. 11). A stellar population observed in this 
phase can have near-IR colors that arc compara- 
ble to those of very old populations. Doubts about 
the validity of the stellar evolution models used in 
the synthesis were raised by OrigUa et al. (1999), 
who were unable to reproduce the colors of RSG 
dominated populous clusters in the Large Magel- 
lanic Cloud with Starburst99. It is therefore of 
interest to evaluate if the use of the Padova tracks 
leads to a significant improvement. 

Since the LMC cluster sample used by Origlia 
et al. (1999) suffers from small-number statistics, 
we decided to use published near-IR photometry of 
blue compact dwarf (BCD) galaxies as test beds. 
BCDs are the archetypal starbursts: they were 
among the galaxies first suggested to experience 
a "burst (or flash) of star formation" by Scarlc, 
Sargent, & Bagnuolo (1973). BCDs are consid- 
ered to be strongly starbursting galaxies powered 
by ionizing stars and an underlying population of 
red stars whose age is still under discussion (e.g., 
Aloisi, Tosi, & Greggio 1999; Thuan, Izotov, & 
Foltz 1999). There is, however, convincing evi- 
dence for a significant number of RSGs in those 
objects with well-established color-magnitude di- 
agrams (e.g., Schulte-Ladbeck et al. 2001). 

A literature search for existing near-IR pho- 
tometry of BCD galaxies resulted in the selection 
of the studies Thuan (1983), Doublier, Caulet, & 
Comte (2001), Noeske et al. (2003), and Telles 
(2004) . These authors obtained J, H, and K large- 
aperture photometry of a statistically significant 
number of BCDs. Thuan's work was done with 
a classical InSb photometer, but through a wide 
aperture of at least 8" , which is large enough to en- 



compass the whole galaxy. The other three studies 
used modern two-dimensional detectors, ensuring 
coverage of the entire underlying galaxies. The 
photometry in all four studies is both sensitive to 
the young starburst, the surrounding field consist- 
ing of potential earlier starburst episodes, and the 
older underlying population. Therefore any RSGs 
present in these galaxies will affect or even domi- 
nate the photometry. 

The photometry values collected from the four 
papers are in Table 5. This table gives the identi- 
fier in col. 1, where KN, VD, TT, and ET stand 
for the initials of the lead authors. Col. 2 lists the 
galaxy names. The near-IR colors [ J — H) and 
{H — K) as well as their corresponding quoted er- 
rors are in cols. 3-6. The table has a total of 
79 galaxy entries, although four are non-unique. 
The photometry was taken as given in the pub- 
lications; no correction for foreground reddening 
was applied by us. This correction was performed 
by Thuan (1983) but not in the other three stud- 
ies. However, the corrections would be less than 
0.03 mag in all cases, which is totally negligible in 
view of other uncertainties. 

The data points from Table 5 are plotted in a 
color-color diagram reproduced in Fig. 20. The er- 
ror crosses are the average of the individual pho- 
tometric errors of each study, off-set in order to 
minimize confusion. Before comparing the data 
to synthetic colors, two corrections need to be ad- 
dressed. First, the J, H and K bands encompass 
several strong nebular emission lines, such as Pa/3, 
which can be strong in the presence of ionizing 
stars (Dale et al. 2004). Vanzi, Hunt, & Thuan 
(2002) obtained near-IR spectra of several galax- 
ies for which they had photometry and determined 
the color corrections that need to be applied to the 
J, H, and K data. Their average corrections arc 
+0.09 and -0.05 mag in the {J - H) and {H - K) 
colors, respectively. These values include both the 
line and the continuous emission. The continuous 
nebular emission (free-free and bound-free) in the 
J, H, and K bands is much more significant, but 
it is accounted for in the synthetic colors and does 
not have to be removed for a comparison. There- 
fore, the pure line emission correction to be ap- 
plied to our sample should be much smaller than 
the above values. We conclude that this effect is 
negligible. 

We also plotted the reddening vector for Ay = 
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0.25 mag in Fig. 20. This corresponds to the av- 
erage internal reddening of BCDs as derived from 
the optical Balmer decrement (Thuan 1983). The 
dispersion, however, can be large. For instance, 
Tol 0610-387 has essentially zero internal redden- 
ing, whereas Mrk 996 has a large value of ~ 0.2 
in {H — K) . The overall effect of both line emis- 
sion and internal reddening is that the bulk of the 
galaxies could move to bluer (smaller) values of 
(J — H) and {H — K) by less than 0.5 mag, with 
a few individual objects moving close to the (0,0) 
point in Fig. 20. 

The synthetic colors predicted for a single stel- 
lar population by Starburst99 and the Geneva 
tracks are plotted in Fig. 20 as the short-dashed, 
long-dashed, and sohd hues for Z = 0.02, 0.004, 
and 0.0004, respectively. The full symbols are tick 
marks to indicate ages of 5, 10, 25, and 100 Myr. 
Before proceeding with a discussion of the impli- 
cations, we introduce Figs. 21, 22, and 23. These 
figures are identical to Fig. 20 in their data con- 
tents but show different model predictions. In 
Fig. 21, we have plotted the predictions from the 
Padova-FTIP-AGB tracks with Z = 0.02, 0.004, 
and 0.0004. Figs. 22 and 23 are the counterparts 
of Figs. 20 and 21, respectively, in that we have 
substituted the predictions for a single starburst 
by those for continuous star formation. 

Figs. 20 and 21 suggest reasonable agreement 
between the bulk of the data and the colors for 
solar chemical composition. The Geneva models 
appear to do marginally better, but this is mainly 
caused by the 10 Myr model point. The predicted 
lifetime of this red loop is quite short so that the 
25 Myr model point should be a better indication 
of the location of the majority of the data. In that 
respect, the Padova and Geneva models are quite 
similar. Of course, a few of the observed BCD will 
in fact have ages of about 10 Myr with little inter- 
nal dust. These objects would require the extreme 
Geneva data point. We terminated the synthetic 
models at an age of 100 Myr. Higher ages are 
unrealistic, as BCDs are defined via their emis- 
sion lines and blue colors. 100 Myr old star bursts 
would not be classified as BCDs anymore. 

The detailed star-formation history of BCDs is 
a subject of debate. Dwarf galaxies are known to 
have had rather complex star-formation histories, 
with periods of quiescence and intermittent bursts 
of star formation (e.g., Greggio et al. 1998). De- 



pending on the burst frequency, the effective star 
formation may be closer to the continuous than 
to the instantaneous case. This is addressed in 
Figs. 22 and 23. As expected, the imprint of the 
RSGs is more diluted, and the colors are more de- 
generate than for a single stellar population. If 
the star-formation histories in BCDs were con- 
stant, the comparison between the data and the 
models would force us to postulate ages far in ex- 
cess of 100 Myr. While this would not necessarily 
be in conflict with observational selections (ion- 
izing stars are continuously replenished), the as- 
sociated gas consiimption would become entirely 
unreasonable. BCDs do not constantly form stars 
over ~ 1 Gyr. Therefore we conclude that the ap- 
propriate star-formation scenario is between those 
plotted in Figs. 20/21 and 22/23, but most likely 
much closer to the instantaneous case of Figs. 20 
and 21. 

Does this suggest consistency between the ob- 
served and synthetic colors? The only track in 
Fig. 20 or 21 that matches the data points is 
the one at solar chemical composition. The other 
tracks at lower abundance are significantly bluer 
and fail to reproduce the observed colors regard- 
less of the assumed age and reddening correc- 
tion. Only the solar models produce RSGs in large 
enough numbers and with sufficiently low T^g to 
intercept with the band delineated by the redden- 
ing corrected data points. A literature survey sug- 
gests that the approximate average oxygen abun- 
dance of our sample is 20% solar. Therefore solar 
composition is entirely excluded, and the applica- 
ble models are those having Z = 0.004. Once we 
are forced to compare the data to the Z = 0.004 
tracks, we arrive at the same conclusion as Origlia 
et al. (1999): The predicted colors of evolutionary 
models for metal-poor populations with a signifi- 
cant RSG component are incorrect. This conclu- 
sion is unchanged for either the Padova or Geneva 
tracks. 

Our residts echo those of Massey & Olsen 
(2003), who demonstrated that both the Geneva 
and Padova evolution models fail to predict the 
location of RSGs and the Large and Small Mag- 
ellanic Clouds. First, the models are too warm at 
any temperature, and second they do not repro- 
duce the smooth decrease in effective temperature 
with increasing luminosity. The higher luminosity 
RSGs are invariably of cooler temperatures, which 
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is just the opposite of the predicted behavior. The 
failure of stellar evolution models, in particular at 
sub-solar chemical composition, implies that any 
spectral synthesis model for metal-poor stellar 
populations with RSGs is doomed to fail as well. 
The evidence of failure at solar chemical composi- 
tion is much weaker, if present at all. 

The disagreement does not entirely come as a 
surprise. Standard stellar models are still based 
on the half-century old mixing-length approxima- 
tion for describing the convective energy trans- 
port. One assumes a fixed value of the ratio of 
mixing length to pressure scale height that directs 
convection in the same manner in all different envi- 
ronments. The mixing-length parameter primar- 
ily governs the stellar radius and is derived from 
a solar calibration. Even if it can be justified in 
Sun- like stars, it may break down in environments 
that are different from that of the Sun. Such en- 
vironments may include not only chemically dif- 
ferent stars but also different evolutionary stages 
of the same star. Although there are efforts to 
find more realistic descriptions of convection (e.g., 
Yi 2003), there is no comprehensive set of stellar 
models based on a more mature convection pre- 
scription. 

Another source of uncertainty is the amount of 
convective overshooting. Overshoot is the pen- 
etrative motion of convective cells, reaching be- 
yond the convective core as defined by the classic 
Schwarzschild criterion. Its effects on stellar evo- 
lution are ubiqitous, most notably on the ratio of 
lifetimes spent in the core hydrogen-burning stage 
and in the shell hydrogen-burning stage. Yi et 
al. (2000) performed an analysis of such effects in 
isochrone and spectrum synthesis. By definition, 
overshooting is effective only for those stars whose 
high masses lead to the formation of a convective 
core. Therefore it affects only young populations, 
producing bluer spectra (Yi 2003). 

Finally, evolution models with rotation are just 
beginning to become available (Maeder & Meynet 
2003). The new models demonstrate the signifi- 
cance of rotation: The additional helium brought 
near the H-burning shell by rotational mixing and 
the larger He-core both lead to a less efficient H- 
burning shell and a smaller associated convective 
zone. Therefore, the stellar radius of rotating stars 
will inflate during the He-burning phase. These 
models account for the formation of more RSGs at 



low Z . with a blue-to- red supergiant ratio in bet- 
ter agreement with observations. It may very well 
be that rotation is a dominant effect and controls 
the interpretation of the spectra of high redshift 
galaxies and the early nucleosynthesis. 

Prustratingly, mixing length, overshooting and 
rotation are the likely culprits for the failure of the 
evolution models, but they cannot be predicted 
from first principles and remain adjustable param- 
eters. Only empirical calibrations can provide con- 
straints. 

Origlia et al. (1999) attempted to derive such 
an empirical calibration. They developed evolu- 
tionary synthesis models from Starburst99 and ap- 
plied them to near-IR spectral features observed 
in the spectra of young Magellanic Cloud clusters. 
The temporal evolution of the first and second 
overtones of CO at 2.29 \vca. (2-0 band head) and 
1.62 \xm (6-3 band head), and of the {U — B), 
{B — V), and (J — K) colors were investigated. 
Origlia et al. utilized star clusters whose ages were 
known independently from spectroscopy, allowing 
them to break the age-related degeneracy. They 
found that the evolutionary models of RSGs with 
sub-solar chemical composition are not reliable for 
any synthesis of the temporal evolution of IR fea- 
tures. While this mirrors the results of Massey & 
Olsen (2003), the integrated light analysis of sin- 
gle clusters allowed them to place constraints on 
the temperature and the fraction of time spent in 
the red part of the (HRD) by massive stars dur- 
ing their core-helium burning phase. The study of 
Origlia et al also addressed the question if a dif- 
ferent choice of the IMF could improve the agree- 
ment, with a negative answer. The near-IR spec- 
tral features from RSGs are practically indepen- 
dent of the IMF. The negligible dependence on 
the IMF simply follows from the fact that RSGs 
provide both the line flux and the adjacent contin- 
uum. Therefore, the IR indices sample a relative 
small stellar mass interval and become IMF inde- 
pendent. 

The key result of Origlia et al. (1999) is that 
the evolution models fail not only in their predic- 
tion of the temperature but also of the lifetimes 
spent in the cool RSG phase. Origlia et al. de- 
rived a set of empirically calibrated spectrophoto- 
metric models by adjusting the RSG parameters so 
that the properties of the observed templates were 
reproduced. Specifically they decreased the effec- 
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tivc temperatures and increased the lifetimes of 
the blue core-helium burning phases. Their study 
suffers from two major limitations: (i) The metal- 
licity range is restricted to that of the Magellanic 
Clouds. Yet, we expect the most dramatic short- 
comings of stellar evolution models at even lower 
metallicity. (ii) An integrated light study is un- 
able to break the mass degeneracy. Origlia et al. 
had to assume that the same correction factors 
apply to RSGs of all masses, which is obviously a 
very drastic assumption with little theoretical or 
empirical justification. 

Therefore we are forced to advice against using 
Starburst99 to model stellar populations at sub- 
solar chemical abundances whose light contains a 
significant RSG component. 

9. Conclusion 

We have enhanced Starburst99 by implement- 
ing a new set of evolutionary tracks which extends 
the useful modeling space of this synthesis code 
to ages older than ~ 1 Gyr and masses below 
0.8 Mq. The prior restrictions in Starburst99 were 
imposed by using the Geneva models, whose are 
optimized for young massive stars. In contrast, 
the newly added Padova tracks are more complete 
and have more detailed input physics at interme- 
diate and low masses. We augmented the Padova 
tracks with dedicated models for AGB stars in or- 
der to complete the evolution of intermediate-mass 
stars until the final thermal pulse prior to entering 
the PN phase. With these enhancements, Star- 
burst99 allows the calculation of synthetic pop- 
ulation properties for virtually any astrophysical 
situation of interest. 

We discussed in this paper the preferred do- 
mains where either the Geneva or the Padova 
models are a better match to observations. The 
former are preferred for young, ionizing popula- 
tions whereas the latter are the models of choice 
when intermediate-age and old population dom- 
inate. The implementation in Starburst99 offers 
users the choice of performing synthesis calcula- 
tions with both sets of tracks in order to help make 
an educated decision as to which are the preferred 
tracks. 

Comparisons with published synthetic proper- 
ties by other groups suggest consistency once dif- 
ferences in the adopted input physics are taken 



into account. In this respect, Starburst99 resem- 
bles GALAXEV most closely, and it agrees very 
well with the predictions of that code. 

We performed a basic sanity check of our mod- 
els by comparing the predicted colors to obser- 
vations of old, intermediate-age, and young clus- 
ters. The agreement is excellent for ages above 
10^ years, which includes the full AGB phase. 
When RSGs dominate, however, models and ob- 
servations at subsolar chemical composition dis- 
agree. This discrepancy is caused by incorrect pa- 
rameters of RSGs in the stellar evolution models. 

All model calculations presented in this paper 

can be repeated by using the on-line version of 
Starburst99 at http://www.stsci.edu/science/starburst99/. 
The upgrades described here are released as ver- 
sion 5.0. In addition, v5.0 has the following 
changes with respect to the currently available 
v4.0. 

Both linear and logarithmic time steps can 
be selected. The latter are a better compromise 
between sufficient time resolution and data volume 
management. 

- We replaced the previously implemented 
Johnson R and I filters with the more widely used 
Cousins R and I filters. Since the Cousins filters 
have shorter effective wavelengths, colors in these 
two passbands tend to be bluer when calculated 
with Starburst99 v5.0. 

- The on-line version of Starburst99 now offers 
more choices for the assumed IMF. While the pre- 
vious version supported a single power law, we im- 
plemented multiple power laws in v5.0. Therefore, 
the commonly used parameterization of Kroupa 
(2002) can be adopted as input. 

We encourage the community to take advantage 
of the improvements whose come with the new re- 
lease of Starburst99. 

We are grateful to Miguel Cerviiio, Cesare 
Chiosi, Daniel Schaerer, and Eva Villaver for pro- 
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this project. Edu Telles kindly made his unpub- 
lished photometric data available to us. This work 
was supported by HST grant GO-9116 from the 
Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Re- 
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Fig. 2.— The Padova tracks with the TIP-AGB 
for (a) Z = 0.02, (b) Z = 0.008, and (c) Z = 0.004. 
The tracks are shifted in log(L/_L0) by 0, 1, 2, ... 
16 dex from their original value of M\ = 0.6, 0.8, 
1.0 .. . 120 Mq, respectively, to improve clarity. 
AGB phases are identified by open triangles. 



Fig. 1. — Geneva (solid) and Padova (dashed) 
tracks for solar composition used in Starburst99 



showing the transition from the intermediate to 
the low-mass regime. The tracks arc shifted in 
logL from their original value by 0, 0.9, 1.6, 
2.1 and 2.7 dex for M = 0.8, 0.9, 1, 1.5, and 
1.7 Mq, respectively. An additional Geneva track 
for 1.7 Mq is plotted to show the transition from 
stars with and without degenerate helium core 
(crosses, shifted by 2.45 from its original value) 
in the Geneva tracks. 
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Fig. 3. — Evolutionary tracks for masses of 1, 3, 
and 5 Mq from Geneva (solid), Padova (dashed), 
and Padova with the TIP-AGB included (dotted), 
(a) Z = 0.02, (b) Z = 0.008, and (c) Z = 0.004. 
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Fig. 4. — Spectral energy distributions for 
young populations with solar composition. Solid: 
Geneva; dashed: original Padova; dotted: 
Padova+TIP-AGB. The far-UV continuum pro- 
duced by the Geneva and Padova models at 3 
Myr differs dramatically and has been labeled 
twice for clarity. Note that the Padova and the 
Padova+TIP-AGB models coincide prior to the 
occurrence of the first AGE stars, producing a 
"dashed-dotted" graph. 




og(Wavelength [Angstroms]) 

Fig. 5. — Spectral energy distributions for 
intermediate-age and old populations. Same as 
for Fig. 4, except for the different age range. The 
spectra for 0.1, 0.5, 1, and 5 Gyr are displaced 
in logL by 2.5, 2.0, 1.5, and 1.0 dex from their 
original value, respectively. 
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Fig. 6. — Evolution of the W-R/0 star ratio. 
Solid: Geneva; dotted: Padova. Solar composi- 
tion. 
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Fig. 7. — Same as Fig. 6, but for the ratio of WC 
over WN stars. 



Fig. 8. — (a) Number of photons with wavelength 
below 504 A at solar composition. Solid line: 
Geneva tracks with the extended atmospheres of 
Smith et al. (2002); long-dashed: Geneva tracks 
with the plane-parallel atmospheres of Lejeune et 
al. (1997; 1998); short-dashed: Padova tracks 
with extended atmospheres; dash-dotted: Padova 
tracks with plane-parallel atmospheres. W-R stars 
first appear at age ~ 3 Myr. (b) Same as in (a) 
but for the number of photons with wavelength 
below 228 A. 



28 




1 -10 

log(Age [Gyr]) log(Age [Gyr]) 




-3 -2-10 1 -3 -3-10 1 

log(Age [Gyr]) log(Age [Cyr|) 



Fig. 9.— MboI, ii500, Mb, My, Mj, and Mk vs. 
time for models with the Geneva (sohd), Padova 
(dashed), and Padova+TIP-AGB (dotted) tracks. 
Luminosity units are magnitudes, except for iisoo- 
The contribution from AGB stars becomes notice- 
able between 0.1 and 3 Gyr at the longest wave- 
lengths. Solar chemical composition. 



Fig. 10.— Contribution of the TIP-AGB phase 
in the Padova tracks to the total light in different 
bands. Solar chemical composition (short dashed 
line), Z = 0.008 (long dashed), and Z = 0.004 
(solid). 
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Fig. 11.— Evolution of {B - V) (top), {V - J) 
(center), and {V — K) (bottom) with time. Left 
panels: Z — 0.02; center panels: Z = 0.008; right 
panels: Z = 0.004. Line types as in Fig. 9. 
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Fig. 13. — Evolution of the total mass loss. Most 
of the mass loss occurs for young populations. 
Models with no mass loss for the TIP-AGB have a 
constant total mass loss for old populations. AGB 
stars provide an incremental increase at very old 
ages. Same notation as before. Solar chemical 
composition. 



log{Age (GyrJ) 

Fig. 12. — Evolution of the mass-loss rate for the 
three sets of tracks. Notice the extension of M 
around log (Age) > —1.3 where AGB stars with 
masses 0.9 < M/Mq < 5.0 have strong mass loss. 
Same notation as before. Solar chemical composi- 
tion. 
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Fig. 14. — The mass-to-light ratio for LboI, Ly, 
and Lk- Same notation as before. Solar chemical 
composition. 
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Fig. 15. — Color evolution for different models 
with solar composition. Bold solid line: this work; 
long-short-dashed: GALAXEV with Padova 
(1994) tracks; long-dash-dotted: GALAXEV with 
Padova (2000) tracks; open small circles: Spectral 
and Padova (1994) isochrones; dotted: Spectral 
and Padova (2000) isochrones; short-dash-dotted: 
Maraston (1998); short-dashed: Fioc & Rocca- 
Volmerange (1997). Starburst99 and GALAXEV 
are plotted together in the left panels because of 
their almost identical physical ingredients. 



Fig. 16. — Near-IR color evolution for solar com- 
position (upper panels) and Z = 0.008 (lower pan- 
els). Bold solid line: this work; long-short-dashed: 
GALAXEV with Padova (1994) tracks; long-dash- 
dotted: GALAXEV with Padova (2000) tracks; 
dotted: Spectral and Padova (2000) isochrones; 
long-dashed: Mouhcine & Langon (2002); short- 
dashed: Fioc & Rocca-Volmerange (1997). 
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Fig. 17. — Evolution ol the mass-to-light ra- 
tio including the full AGB phase for Lbo\, Ly 
and Lk at solar composition. Upper set of pan- 
els: fixed mass for all models except Maraston 
(1998); lower set: variable mass for all models 
except Fioc & Rocca-Volmerange (1997). Solid 
line: this work; long-short-dashed: GALAXEV 
with Padova (1994) tracks; long-dash-dotted: 
GALAXEV with Padova (2000) tracks; short- 
dash-dotted: Maraston; short-dashed: Fioc & 
Rocca-Volmerange. 



Fig. 18.— {V - I) vs. {B - V) two-color di- 
agram of NGC 5128. The colors predicted by 
the models are indicated by the three curves for 
Z = 0.0004 (solid), 0.004 (long-dashed), and 0.02 
(short-dashed). Tick marks indicate ages of 0.1, 
1, 5, and 10 Gyr. The points (together with the 
schematic error cross at 0,1.5) represent the globu- 
lar cluster photometry of Peng et al. (2004a). The 
arrow is the reddening vector for a foreground red- 
dening of E{B -V) = 0.115. 
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Fig. 19. — Same as Fig. 18, but for the star clus- 
ters in NGC 4038/39 of Whitmore et al. (1999). 
Open stars: old globular clusters; open triangles: 
intermediate-age clusters; open squares: young 
clusters. No reddening correction was applied. 
The models are identical to those in Fig. 18. 
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Fig. 20. — Color-color diagram of four samples 
of BCDs compared with starburst models of Z = 
0.001 (sohd hue, blue), 0.004 (long-dashed line, 
green), and solar Z = 0.02 (short-dashed line, red) 
with Geneva tracks. Filled symbols indicate the 
ages of the models. Open symbols show the data: 
triangles (Doublier et al.; 2001), squares (Noeske 
et al.; 2003), starred circles (Thuan; 1983), and 
stars (Telles; 2004). The error bars are the dis- 
persion value in each sample. The solid vector 
indicates the reddening correction for Ay = 0.25 
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Fig. 21. — Same as Fig. 20 for starburst models of 
Z = 0.0004 (solid line, blue), 0.004 (long-dashed 
line, green), and Z — 0.02 (short-dashed hne, red) 
with Padova tracks. 
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Fig. 23. — Same as Fig. 20 for continuous star for- 
mation models of Z = 0.0004 (solid line, blue), 
0.004 (long-dashed hne, green), and Z = 0.02 
(short-dashed line, red) with Padova tracks. 



Fig. 22. — Same as Fig. 20 for continuous star 
formation models oi Z — 0.001 (solid line, blue), 
0.004 (long-dashed hne, green), and Z = 0.02 
(short-dashed line, red) with Geneva tracks. 
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Table 1 

Adopted solar abundances of the main elements in the Geneva and Padova tracks. 



Atomic Number Element Mass Fraction 



1 


H 


7.06 X 


10-1 


2 


He 


2.75 X 


10-1 


3 


Li 


9, 


,35 X 


10-3 


4 


Be 


1.66 X 


10-1" 


5 


B 


4, 


,73 


X 


10-9 


6 


C 


3, 


,03 


X 


10-^ 


7 


N 


1, 


,11 


X 


10-^ 


8 


O 


9, 


,59 


X 


10-^ 


9 


F 


4, 


,05 


X 


lo-'^ 


10 


Ne 


1, 


,62 


X 


10-^ 


11 


Na 


3, 


,34 


X 


10-^5 


12 


Mg 


5 


,15 


X 


10-* 


13 


Al 


5 


,81 


X 


10-"^ 


14 


Si 


6, 


,53 


X 


10-* 


16 


S 


3, 


,96 


X 


10-* 


25 


Mn 


1 


,33 


X 


10-^ 


26 


Fe 


1, 


,17 


X 


10-^ 
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Table 2 

Mass fractions of H, He, C, N, and O used in the Geneva and Padova tracks. 







Padova Tracks 






Element 


Z = 0.0004 


Z = 0.004 


Z = 0.008 


Z = 0.02 


Z = 0.05 




7.70(-01) 


7.56(-01) 


7.42(-01) 


7.00(-01) 


5.98(-01) 




2.30(-01) 


2.40(-01) 


2.50(-01) 


2.80(-01) 


3.52(-01) 




9.90(-05) 


9.89(-04) 


1.98(-03) 


4.94(-03) 


1.24(-02) 




2.50(-05) 


2.48(-04) 


4.97(-04) 


1.24(-03) 


3.10(-03) 




2. 11 (-04) 


2.11(-03) 


4.23(-03) 


1.06(-02) 


2.64(-02) 






Geneva Tracks 








Z = 0.001 


Z = 0.004 


Z = 0.008 


Z = 0.02 


Z = 0.04 




7.56(-01) 


7.44(-01) 


7.28(-01) 


6.80(-01) 


6.20(-01) 


^He 


2.43(-01) 


2.52(-01) 


2.64(-01) 


3.00(-01) 


3.40(-01) 


= C 


2.43(-04) 


9.73(-04) 


1.95(-03) 


4.87(-03) 


9.73(-03) 


''N 


6.20(-05) 


2.47(-04) 


4.95(-04) 


1.24(-03) 


2.47(-03) 


«o 


5.27(-04) 


2.11(-03) 


4.22(-03) 


1.05(-02) 


2.11(-02) 
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Table 3: Properties of the two sets of tracks. 



Property 


Geneva 


Padova+TIP-AGB 




0.05, 0.2, 0.4, 1, 2 


0.02, 0.2, 0.4, 1, 2.5 


Upper Mass 


120 


120 


Lower Mass 


0.8 


0.15 


Grid Poiiils 


51 


84 


PllilSCS 


AlH, IKili, ZAllB" , and TP-AC;b'' 


MS, H.GB. IIB, TP-AGB, and Tli>-AGB' 



"The HB phase for intermediate-mass stars in Geneva models included only at the beginning of the helium burning on the 

zero-age horizontal branch (ZAHB). 
''For stars with M < 1 Mq, the tracks have the RGB phase for Z = 0.004 and 0.008. Otherwise the evolution is limited to the 

early evolution after the MS. The RGB is the last jjliase modeled in stars with 1 < M/Mq < 1.7. 
'^For stars with 0.6 < M/Mq < 0.9, the tracks follow the evolution up to the RGB phase only. 
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Table 4 

AGB IMPLEMENTATION IN StARBURST99 AND OTHER EVOLUTIONARY SYNTHESIS CODES. 



Designation 


Stellar models 


AGB phase 


Mm ax 




SEDs 


Photometry 


Reference 


Starburst99 


Padova (2000) tracks 


Vassiliadis & Wood (1993) 


5 


Leje 


une et al. (1997;1998) 


Johnson-Cousins 


This work 


GALAXEV (2003) 


Padova (2000) tracks 


Vassiliadis & Wood (1993) 




Le 


Borgne et al. (2003) 


Johnson-Cousins 


Bruzual & Chariot (2003) 


GALAXEV (2003) 


Padova (1994) tracks 


Vassiliadis & Wood (1993) 




Le 


Borgne et al. (2003) 


Johnson-Cousins 


Bruzual & Chariot (2003) 


Spectral 


Padova (2000) isochrones 


Groenewegcn & de Jong (1993)^ 




Leje 


jine et al. (1997; 1998) 


Johnson-Cousins 


Vazquez et al. (2003) 


Spectral 


Padova (1994) isochrones 


Boothroyd & Saekmann (1988)'=' 




Leje 


jne et al. (1997; 1998) 


Johnson- Cousins 


Vazquez et al. (2003) 


Mouhcine & Langon 


Padova (1994) tracks 


Wagenhuber & Groencwegen (1998)'= 


7 


La 


ngon & Wood (2000) 


Johnson- Cousins 


Mouhcine & Langon (2002) 


Maraston 


Castellani et al. (1992) 


Blocker & Schonberner (1991) 


8 




Kurucz (1979) 


Johnson 


Maraston (1998) 


PEGASE 


Padova (1994) tracks 


Schonberner (1983), Blocker (1995) 


6 


Gu 


nn & Stryker (1983) 


Johnson-Cousins 


Fioc Rocca-Volmerange (1997) 



^The mass loss prescription of Vassiliadis &c Wood (1993) is used. 
^The Reimers (1975) formula for the mass-loss rate is used. 
^The Blocker (1995) formula for the mass loss rate is used. 



Table 5 

Sample of BCD galaxies used for comparison with our models. 



Id 


Object 


(J - H) 


A(J - H) 


(H - K) 


A(H - K) 


IKN 


Tol 3 


0.60 


0.07 


0.23 


0.06 


2KN 


Tol 1400-41 Ic 


0.36 


0.08 


0.24 


0.09 


3KN 


Pox 4c 


0.39 


0.13 


0.29 


0.13 


4KN 


Pox 4B 


0.40 


0.15 


0.20 


0.15 


5KN 


UM 448 


0.65 


0.10 


0.20 


0.11 


6KN 


UM 461 


0.47 


0.06 


0.20 


0.06 


7KN 


He 2-10 


0.57 


0.06 


0.24 


0.07 


8KN 


IC 4662c 


0.59 


0.12 


0.12 


0.12 


9KN 


Mrk 178 


0.45 


0.15 


0.24 


0.16 


lOKX 


Mrk l.'i2n 


O.fiCi 


0.09 


n.ni 


n.in 


IIKA 


lla.ro 14 


(1.71 


(1.12 


0.1(1 


0.12 


12VD 


Haro 14 


0.60 


0.24 


0.18 


0.14 


13VD 


Mrk 996 


0.47 


0.24 


0.42 


0.14 


14VD 


Tol 0610-387 


0.44 


0.24 


0.18 


0.14 


15VD 


Tol 0957-279 


1.02 


0.24 


0.03 


0.14 


16VD 


Fairall 301 


0.56 


0.24 


0.29 


0.14 


17VD 


Tol 0954-293 


1.00 


0.24 


-0.02 


0.14 


18TT 


Mrk 36 


0.38 


0.14 


0.42 


0.15 


19TT 


Mrk 49 


0.54 


0.14 


0.25 


0.15 


20TT 


Mrk 59 


0.40 


0.14 


0.44 


0.15 


21TT 


Mrk 86 


0.69 


0.03 


0.20 


0.04 


22TT 


Mrk 89 


0.72 


0.06 


0.07 


0.07 


23TT 


Mrk 116 


0.52 


0.18 


0.11 


0.27 


24TT 


Mrk 140 


0.70 


0.04 


0.24 


0.06 


25TT 


Mrk 169 


0.67 


0.03 


0.35 


0.04 


26TT 


Mrk 209 


0.40 


0.08 


0.10 


0.12 


27TT 


Mrk 297 


0.75 


0.04 


0.29 


0.06 


28TT 


Mrk 313 


0.80 


0.04 


0.26 


0.06 


29TT 


Mrk 314 


0.50 


0.08 


0.09 


0.08 


30TT 


Mrk 324 


0.63 


0.05 


0.19 


0.07 


31TT 




0.64 


0.06 


0.21 


0.07 


32TT 


Mrk 370 


0.61 


0.04 


0.26 


0.06 


33TT 


Mrk 401 


0.75 


0.03 


0.26 


0.04 


34TT 


Mrk 527 


0.83 


0.04 


0.31 


0.06 


35TT 


Mrk 600 


0.54 


0.11 


0.07 


0.11 


36TT 


I-Zw 123 


0.54 


0.15 


0.24 


0.15 


37TT 


II-Zw 70 


0.55 


0.04 


0.17 


0.05 


38TT 


II-Zw 71 


0.61 


0.04 


0.15 


0.06 


39TT 


III-Zw 102 


0.80 


0.04 


0.29 


0.06 


40TT 


Haro 2 


-0.01 


0.07 


0.25 


0.04 


41TT 


Haro 3 


-0.08 


0.04 


0.17 


0.04 


42TT 


Haro 9 


0.59 


0.14 


0.14 


0.14 


43TT 


Haro 14 


0.61 


0.05 


0.11 


0.06 


44TT 


Haro 15 


0.59 


0.04 


0.21 


0.06 


45TT 


Haro 18 


0.63 


0.05 


0.10 


0.07 


46TT 


Haro 20 


0.64 


0.07 


0.22 


0.07 


47TT 


Haro 21 


0.66 


0.05 


0.16 


0.06 


48TT 


Haro 23 


0.56 


0.04 


0.18 


0.08 


49TT 


Haro 28 


0.67 


0.14 


0.18 


0.14 


50TT 


Haro 36 


0.55 


0.15 


0.08 


0.17 


51TT 


Haro 38 


0.55 


0.15 


0.14 


0.15 


52TT 


Haro 43 


0.71 


0.30 


0.39 


0.29 


53ET 


C0828a 


0.73 


0.05 


0.35 


0.05 


54ET 


II-Zw 40 


0.76 


0.05 


0.14 


0.05 


55ET 


I-Zw 18 


0.24 


0.05 


0.33 


0.05 


56ET 


Mrk 36 


0.22 


0.05 


0.31 


0.05 


57ET 


UM 448 


0.63 


0.10 


0.30 


0.10 


58ET 


UM 461 


0.03 


0.10 


0.01 


0.10 
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Table 5 — Continued 



Id Object (J - H) 


A{J -H) 


(H-K) 


A{H - K) 


59ET UM 559 0.14 


0.10 


-0.05 


0.10 



40 



